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CHAPTER  I.  INTRODUCTION 


The  scientific  and  technological  developments  in  all  fields  tend  toward 
equipment  of  ever-increasing  complexity.  This  complexity  brings  with  it  the 
inescapable  concomitant  —  delicacy.  In  addition,  the  multinational  posture  of 
many  of  the  leading  corporations  demand  that  this  delicate  equipment  be 
transported  to  virtually  every  corner  of  the  globe.  Consequently,  extremely 
sophisticated  container  systems  are  required  to  protect  this  delicate  equip¬ 
ment  during  its  journey  throughout  the  world. 

Many  items  come  adequately  packaged  by  nature.  Consider  the  grape¬ 
fruit  with  its  outside  protective  shield  which  protects  it  from  damage  on  the 
tree  and  in  transit  to  the  consumer’s  table.  Unfortunately,  many  of  man’s 
physical  creations  do  not  possess  a  protective  shield,  since  to  provide  one 
would  be  economically  infeasible.  If  a  protective  shield  capable  of  preventing 
any  damage  to  a  missile  system  were  incorporated  into  the  basic  structure  of 
the. missile,  the  missile  would  be  incapable  of  flight  due  to  the  weight  burden. 
Consequently,  missile  systems  are  designed  with  only  an  in-flight  protective 
shield  incorporated.  The  container  designer  is  then  assigned  the  task  of 
designing  a  protective  shield,  for  the  various  segments  of  the  missile  system, 
to  be  utilized  for  ground  handling  protection.  The  fundamental  product  of  the 
design  effort  is  identified  as  a  container  system.  The  container  system  is  the 
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dev  utilized  to  protect  the  item  during  transit  to  its  destination  point, 
hopefully  without  damage  and  at  a  minimum  cost. 

In  an  attempt  to  satisfy  the  demand  for  unique  container  systems,  the 
container  designer  is  confronted  with  a  three  faceted  problem.  He  is  expected 
to  protect  items  of  greater  and  greater  fragility  against  greater  and  greater 
environmental  shock  hazards,  at  lower  and  lower  values  of  weight,  volume, 
and  cost.  In  response,  the  container  designer  has  turned  to  a  variety  of  bulk 
cushioning  materials,  primarily  because  of  the  rapidly  mounting  costs  for  the 
design  of  special  shock-isolation  systems.  However,  to  provide  protection 
against  the  anticipated  hazards  encountered  in  transit,  the  container  designer 
requires  information  on  the  physical  properties  of  cushioning  materials,  and 
the  methods  of  configurations  available  for  the  most  economical  use  of  the 
cushioning  materials. 

The  Container  Cushioning  Concept 

.Cushioning  systems  are  incorporated  into  container  systems  to  protect 
fragile  items0  The  degree  of  fragility  of  the  item  determines  the  amount  of 
cushioning  required  to  protect  it  from  damage  during  handling  and  shipment. 
Some  items  are  inherently  strong  and  rugged  except  for  one  or  more  exceed¬ 
ingly  fragile  components.  When  the  fragile  components  cannot  be  removed 
for  separate  cushioning,  the  entire  item  must  be  treated  as  fragile,  even 
though  the  result  may  be  an  unavoidably  large,  cumbersome  container. 

A  basic  container  cushioning  problem  is  shown  in  Figure  1  wh  ch 
utilizes  a  foamed  material  as  the  cushioning  system.  The  item  which  is  to  be 
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Figure  1.  Container  cushioning  system. 

protected  from  damage  is  located  in  the  center.  A  rigid,  wooden  outside 
container  encases  the  protected  item  and  the  foamed  cushion.  The  thickness 
of  the  cushioning  material  is  apparent.  In  this  simplified  sketch,  when  the 
container  is  accidentally  dropped,  it  falls  freely  until  a  rigid  surface  of  some 
type  is  encountered.  Assuming  a  nonresilient  rigid  surface  which  the  con¬ 
tainer  strikes  in  a  flat  fashion,  it  is  noted  that  the  container,  being  rigid, 
impacts  with  considerable  force.  However,  if  the  proper  foamed  cushion  has 
been  provided,  the  item  to  be  protected  does  not  stop  as  quickly  as  the  outside 
container  (from  a  tim<~  point  of  view)  since  the  protected  item  compresses 
into  the  foamed  cushioning  material.  If  the  cushion  is  thick  enough  and  of 
proper  density  for  the  ensuing  impact,  the  protected  item  will  not  be  damaged. 
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The  challenge  then  is  to  determine  the  proper  size  and  density  of  a  foamed 
cushioning  material  to  utilize  as  a  shock  mitigation  system. 

The  container  system  includes  the  container  and  the  cushioning  system, 
regardless  of  the  cushioning  material  utilized.  The  cushioning  system  relates 
to  the  cushioning  material  itself,  together  with  the  cushioning  configuration 
utilized.  In  Figure  1  the  cushioning  configuration  shown  utilizes  side  pads. 
Since  various  outer  containers  can  be  devised  to  encase  almost  any  size  of 
item,  the  crucial  determination  reduces  to  cushion  thickness  and  cushion 
configuration. 


Research  Objective 

The  objective  of  this  research  effort  is  to  develop  the  methodology  for 
confined,  temperature  sensitive,  cushioning  system  modeling.  These  models 
should  be  capable  of  predicting  the  dynamic  cushioning  performance  of  the 
Minicel  cushioning  material  in  the  confined  state.  A  secondary  objective  is  to 
assess  the  relationship  between  unconfined  drop  test  results  on  the  Minicel 
material. 

In  response  to  this  research  objective,  a  systematic  study  of  back¬ 
ground  material  was  conducted.  The  evolution  of  cushion  design  is  discussed 
in  Chapter  n.  The  experimental  design  necessary  for  drop  test  data  acquisi¬ 
tion  is  given  in  Chapter  III.  The  parameters  of  the  overall  impact  response 
model  are  identified,  together  with  the  viscoelastic  theory  upon  which  the 
model  is  predicated,  in  Chapter  IV. 


The  development  of  the  interior  box  model,  together  with  the  total  box 
model,  is  presented  in  Chapter  V.  This  includes  the  individual  dynamic 
cushioning  equations  for  the  interior  box,  the  general  model  for  the  interior 
box,  the  individual  total  box  equations,  and  the  general  model  for  the  total 
box.  Chapter  VI  contains  the  results  of  the  validation  of  the  two  general 
models.  Chapter  VII  combines  the  two  models  for  an  analysis  of  the  tempera¬ 
ture  effects  upon  the  cushioning  material.  Chapter  VIII  contains  the  conclu¬ 
sions  and  recommendations  resulting  from  the  research  effort. 

It  is  anticipated  that  the  development  of  the  overall  container  impact 
response  mathematical  model  will  result  in  a  more  realistic  predictor  of 
cushioning  response  within  the  container  and,  in  turn,  will  prove  to  be  more 
cost  effective  to  cushioning  designers. 


CHAPTER  H.  THE  EVOLUTION  OF  CUSHION  DESIGN 


Prior  to  the  1940’s,  proper  cushioning  was  determined  via  the  ’’cut  ard 
try”  method.  One  used  his  judgment  in  selecting  the  optimum  cushioning 
material,  the  cushion  configuration,  and  the  material  thickness.  If  it  pro¬ 
vided  the  required  protection  it  was  deemed  adequate. 

Early  scientific  investigation  took  the  classic  stress-strain  approach. 
By  slowly  applying  a  load  to  a  specimen  and  measuring  the  resulting  deflection, 
data  were  collected  for  the  plotting  of  strain  curves  for  various  cushioning 
materials.  These  curves  provided  information  pertaining  to  the  amount  of 
energy  stored  by  the  specimen  and  the  point  at  which  resilient  properties  were 
lost  by  the  specimen. 

Prompted  by  the  increased  sophistication  of  equipment  and  the  advances 
in  cushioning  material  fabrication,  serious  scientific  research  into  cushioning 
design  was  begun.  R.  D.  Mindlin  of  Columbia  University  is  credited  with 
having  made  the  first  scientific  effort  in  this  field  in  1945. 

Mindlin  [1]  first  derived  equations  of  motion  by  treating  the  cushion  as 
a  linear  spring.  Then,  assuming  a  linear  load-displacement,  he  used  the 
equations  of  motion  to  find  the  maximum  displacement  of  the  packaged  item 
within  the  container  to  insure  the  item  does  not  attain  a  critical  acceleration. 
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Since  in  practice  a  container  system  rarely  has  linear  spring  characteristics, 
Mindlin  developed  equations  to  describe  the  maximum  displacement  for  non¬ 
linear  springs. 

In  determining  the  optimum  cushion  thickness,  Mindlin  used  the  maxi¬ 
mum  allowable  acceleration,  the  predicted  drop  height,  and  the  weight  of  the 
item  to  find  the  necessary  displacement  of  the  cushion.  The  displacement 
distance  can  be  compared  with  the  load-displacement  data  for  the  various 
cushioning  materials  from  which  a  selection  can  be  made. 

In  this  approach,  load-displacement  data  were  acquired  by  applying 
successively  increasing  loads,  with  weights  or  in  a  load  testing  machine,  and 
measuring  the  corresponding  displacements.  This  type  of  testing  is  cate¬ 
gorized  as  static  because  of  the  low  loading  rate. 

Mindlin’s  development  of  the  mathematical  equations  describing  the 
maximum  acceleration  and  displacement  of  the  protected  item  is  considered 
classical  and  his  approach  to  cushion  design  was  not  significantly  changed  for 
over  a  decade. 

Although  information  on  the  static  properties  of  cushioning  materials 
is  helpful,  it  does  not  permit  the  cushioning  designer  to  predict  the  expected 
limit  of  protection  for  an  item  under  severe  handling  conditions.  Protective 
cushioning  must  be  capable  of  responding  to  a  sharp  force  of  short  duration  if 
adequate  item  protection  is  to  be  provided. 

Figure  2,  taken  from  Mindlin’s  paper,  depicts  an  idealized  mechanical 
system  representing  an  item  during  a  drop  test.  The  outer  container  is 


Figure  2.  Idealized  mechanical  system  representing  item 
during  a  drop  test  [1]. 


represented  by  m^,  and  m 2  is  the  protected  item.  The  cushion  between  the 
item  and  the  outer  container  is  represented  by  a  spring,  . 

In  addition  to  Mindlin,  extensive  research  effort  was  conducted  by 
K.  Q.  Kellicutt,  R.  E.  Jones,  and  D.  L.  Hunzicker  at  the  Forest  Products 
Laboratory  of  the  U.  Department  of  Agriculture.  Another  active  researcher 
was  A.  M.  Underhill  oi  the  General  Electric  Company.  Most  of  these  efforts 
were  directed  toward  collection  of  load-displacement  data  on  the  various 
materials  available  at  that  time  for  application  to  the  cushion  design  methods 
of  Mindlin  [1]  and  Masel  [2], 

Gretz  [3]  extended  the  force-displacement  approach  to  derive  an 
efficiency  factor  for  the  cushioning  material.  The  efficiency  factor  was 
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then  used  in  the  following  equation  to  arrive  at  the  necessary  cushion  thick¬ 
ness: 


Cushion  Thickness  = 


(Drop  Height)  ( Cushion  Deflection)  (Efficiency) 
(Maximum  G’s) 


Janssen  [4]  utilized  the  applied  load  to  the  cushion  and  the  energy 
absorbing  capacity  of  the  cushion  (defined  as  the  area  under  the  load- 
displacement  curve)  to  form  a  dimensionless  ratio  which  he  called  the  cushion 
factor  or  "J”  factor.  The  cushion  factor  was  the  ratio  of  optimum  stress  to 
optimum  strain  and  was  used  to  calculate  the  cushion  thickness. 

These  cushion  design  methods  were  adequate  for  most  of  the  cushioning 
materials  available  during  the  1940’s  and  early  1950’s.  However,  with  the 
development  and  proliferation  of  a  myriad  of  foam  materials  in  the  late  1950’s, 
it  became  apparent  to  cushion  designers  that  the  load-displacement  curves 
produced  through  static  testing  did  not  reflect  the  true  nature  of  viscoelastic 
materials. 

Consequently,  it  was  necessary  to  perform  dynamic  tests  to  investigate 
foamed  material  performance.  A  dynamic  test  is  performed  by  having  a 
plate,  whose  weight  may  be  varied  to  achieve  different  static  stresses,  drop 
onto  the  cushioning  material.  The  plate  is  instrumented  to  sense  the  peak 
acceleration  it  attains.  Figure  3  shows  the  apparatus  used  in  this  type  of  test. 
This  type  of  testing  is  referred  to  as  flat  pad  or  unconfined  dynamic  testing. 

Kerstner  [  5]  pointed  out  the  inconsistencies  found  in  the  dynamic  test¬ 
ing  of  only  one  thickness  of  cushion  and  then  utilizing  that  data  to  predict  the 
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behavior  of  other  thicknesses.  He  also  stated  that  the  cushion  factors  would 
vary  according  to  the  area  of  the  cushion  specimen  tested.  He  attributed  this 
to  materials  which  had  a  relatively  high  pneumatic  characteristic. 

Stern  (6J  described  the  phenomenon  of  the  drastic  difference  between 
the  response  of  some  materials  to  dynamic  loading  and  that  of  static  loading. 

He  attributed  this  to  the  pneumatic  effect  of  the  partially  entrapped  air  in  the 
material.  The  effect  was  found  to  be  most  pronounced  in  interconnected 
cellular  materials  and  varied  with  the  dimensions  of  the  material.  In  a  static 
test  the  air  had  sufficient  time  to  escape  without  absorbing  energy.  In  a  dynamic 
test  the  load  approached  that  of  a  free  falling  object  and  the  rate  of  loading 
upon  the  cushion  was  much  more  rapid.  As  the  cushion  was  compressed,  the 
entrapped  air  in  the  cellular  material  was  not  able  to  escape  as  the  cells 
suddenly  collapsed. 

Figure  4  shows  the  difference  in  energy  absorption  ability  during  static 

and  impact  loading.  E  is  the  energy  dissipated,  E  is  the  energy  stored  in 

fl  s 

the  cushion,  and  e  is  the  maximum  strain  the  cushion  encounters  [71.  The 
m 

open  circles  represent  static  loading  while  the  solid  circles  depict  the  ratio  of 
dissipated  energy  to  stored  energy  of  a  cushion  under  impact  loading  conditions. 

Kerstner’s  solution  to  the  dilemma  of  static  versus  dynamic  data  was 
to  use  peak  acceleration  versus  static  stress  graphs.  The  peak  accelerations 
obtained  during  dynamic  testing  were  plotted  against  the  change  in  load  or 
static  stress.  This  resulted  in  the  graphs  shown  in  Figure  5  [8],  Peak 
acceleration-static  stress  graphs  were  developed  for  various  cushioning 
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Figure  4.  The  ratio  of  dissipated  energy  to  stored 
energy  during  static  and  impact  loading  [7], 

materials  at  various  drop  heights.  The  fragility  level  of  the  item  to  be  pro¬ 
tected  was  determined  by  anticipating  the  maximum  acceleration  the  item 
could  undergo  without  being  damaged  or  destroyed.  The  graphs  are  read  by 
finding  the  fragility  level  on  the  vertical  axis  and  drawing  a  horizontal  line 
from  that  point,  intersecting  the  dynamic  cushioning  curves.  The  cushion 
thickness  associated  with  the  first  curve  which  falls  below  the  line,  within  the 
expected  stress  levels  the  item  will  encounter,  is  the  proper  cushion  amount 


to  utilize 
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Although  Kerstner  indicated  that  his  method  could  also  be  utilized  to 
produce  the  cushion  efficiency  factor,  the  dynamic  cushioning  curves  gained 
lone  acceptance  and  are  still  utilized  today. 

Extensive  research  work  has  been  underway  over  the  past  two  decades 
resulting  in  the  development  of  dynamic  cushioning  curves  for  various  cush¬ 
ioning  materials.  However,  many  of  these  early  curves  were  derived  from 
data  collected  at  ambient  temperature. 

McDaniel  and  Wyskida  [  9  J  were  the  first  to  conduct  extensive  research 
on  the  effects  of  temperature  upon  cushioning  materials.  It  was  demonstrated 
conclusively  that  the  temperature  effects  were  significant.  In  1975,  McDaniel 
[  10  ]  derived  the  first  general  mathematical  model  to  predict  peak  acceleration 
for  flat  pad  drops  at  different  temperatures.  He  generated,  via  drop  tests, 
dynamic  test  data  on  Minicel,  the  commercial  name  for  a  closed  cell,  cross- 
linked  polyethylene  foam  material.  Using  a  stepwise  regression  procedure 
he  was  able  to  develop  the  general  mathematical  model  from  the  experimental 
data.  McDaniel’s  general  model  has  been  used  to  develop  specific  mathe¬ 
matical  models  for  polyester  and  polyether  types  of  polyurethane,  and  Dow 
Ethafoam  polyethylene. 

F  igure  6  is  a  typical  set  of  dynamic  cushioning  curves  used  by  modern 
package  designers.  The  three  curves  represent  three  different  temperatures. 
Points  on  the  curves  correspond  to  the  peak  acceleration  (given  in  G’s) 
experienced  by  the  cusnioned  item  at  the  corresponding  static  stresses.  The 
horizontal  line  depicts  the  fragility  level  to  which  the  item  under  consideration 
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Figure  6.  Typical  dynamic  cushioning  curves  at  selected 
temperatures  derived  from  unconfined  test  data. 

must  be  protected.  The  package  designer  must  choose  a  material  with  a 
material  thickness  which  possesses  dynamic  cushioning  curves  for  the  tem¬ 
perature  range  for  which  he  is  designing  and,  furthermore,  a  portion  of  the 
curves  must  be  below  the  specified  fragility  level.  A  separate  set  of  curves 
is  necessary  for  each  unique  set  of  parameters  (i.e. ,  drop  height,  cushion 
thickness,  and  material) . 

Wyskida,  Johannes,  and  Wilhelm  [11]  also  devised  a  cushion  optimiza¬ 
tion  program  on  the  Hewlett-Packard  (HP)  981 5A  with  a  plotting  capability. 
The  modern  design  engineer  has  access  to  this  program  which  includes  five 


cushioning  materials 
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As  stated  previously,  the  dynamic  cushioning  curves  used  by  cushion 
designers  are  derived  from  data  obtained  when  drop  tests  are  made  on  a  flat 
pad  (unconfined  testing) .  Even  the  impact  response  model  developed  by 
McDaniel  [10]  was  derived  from  flat  pad  drops.  The  cushion  itself  is  not 
enclosed  in  any  way.  In  spite  of  this  advancement  of  the  state-of-the-art, 
many  items  must  be  placed  within  a  container  for  shipment  purposes. 

McDaniel  did  not  consider  the  effects  of  a  container  in  his  model  development. 
This  lack  of  data  on  container  effects  and  inability  to  predict  them  has  been  a 
major  concern  to  cushion  designers. 

Since  the  flat  pad  impact  response  has  been  modeled,  it  is  now  possible 
to  extend  it  to  the  response  within  a  container.  Although  there  have  been  very 
few  direct  comparisons  in  this  area,  the  assumption  has  been  that  the  flat -pad 
data  result  in  conservative  designs.  Grabowski  [12],  using  foamed  poly¬ 
styrene  and  flat  drops  in  a  corrugated  board  container,  wood  crate,  and  steel 
cylinder  at  30-in.  drop  heights,  found  that  the  peak  accelerations  were  lower 
than  those  for  corresponding  static  stresses  and  thicknesses  of  the  flat  pad 
test.  If  the  impact  response  of  the  overall  container  system  (i.e. ,  the  con¬ 
tainer  together  with  the  cushioning  material  as  one  unit)  could  be  modeled,  it 
would  undoubtedly  result  in  substantial  savings  to  the  cushion  design. 

Prior  to  modeling  the  overall  container  system,  it  is  necessary  to 
develop  an  experimental  design  which  is  capable  of  acquiring  the  experimental 
data  for  modeling  purposes.  This  is  the  topic  of  Chapter  III. 


CHAPTER  III.  EXPERIMENTAL  CONSIDERATIONS 

Dynamic  cushioning  curves  are  derived  from  data  obtained  when  drop 
tests  are  conducted  on  a  flat  pad  cushion.  This  method  of  obtaining  impact 
response  Was  not  designed  to  take  into  consideration  the  effects  of  an  outside 
container  upon  the  impact  response.  Mazzei  [13]  found  that  there  was  a 
definite  difference  in  confined  (outside  container)  and  unconfined  (no  container) 
test  results.  He  attributed  this  difference  to  pneumatic  effects  within  the 
container.  Since  there  are  known  to  be  container  effects,  the  lack  of  data  on 
container  effects  and  the  inability  to  predict  them  has  been  a  major  concern  to 
cushion  system  designers.  Consequently,  the  immediate  objective  is  tc  deter¬ 
mine  whether  the  methodology,  developed  by  McDaniel  in  constructing  a  mathe¬ 
matical  model  for  impact  response  from  unconfined  data,  could  be  extended  to 
the  development  of  a  general  model  for  data  collected  from  drop  tests  of  an 
overall  container.  Consequently,  a  drop  test  experiment  was  designed  from 
which  data  were  available  to  apply  these  mathematical  modeling  techniques. 

McDaniel  [10]  developed  a  specific  model  for  the  impact  response  of 
the  Hercules  Minicel  cushioning  material  based  upon  flat  pad  cushion  data. 
Minicel  possesses  the  ability  to  withstand  extreme  temperatures,  without  a 
degradation  in  cushioning  ability.  Consequently,  a  2-in.  thick  Minicel 
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cushion  was  selected  as  the  bulk  cushioning  material  for  use  in  the  overall 
container  experimental  model. 

Physical  Description 

Twenty-five  complete  container  systems  were  prepared  for  experi¬ 
mental  testing.  A  complete  container  system  is  shown  in  Figure  7.  The  two 
major  components  are  the  outside  container  and  the  interior  box.  The  interior 
box  is  a  19-in.  plywood  cube  protected  by  a  2-in.  thickness  of  Minicel  cushion¬ 
ing  material  configured  as  corner  void  pads.  This  interior  box  was  enclosed 
in  a  military  standard  cleated  plywood  shipping  container  as  shown  in 
Figure  7. 


Figure  7.  interior  box  and  outside  container. 
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Six  different  static  stress  levels  were  obtained  by  varying  the  total 
area  of  the  cushioning  material  on  the  six  faces  of  the  interior  box  so  that 
each  face  of  the  box  represents  a  different  stress  level.  The  relationship 
between  static  stress  a,  weight  W  ,  and  total  cushion  area  A  ,  was  used  to 
determine  the  size  of  the  corner  pads  in  the  following  equation: 

W 


Using  the  average  weight  of  the  interior  boxes,  25.5  lb,  and  the  following 


cushion  pad  dimensions,  the  resulting  static  stress  levels  are: 


Cushion  Pad 

Total 

Static  Stress 

Side 

Area 

4. 

Levels 

(in.)  (in.2)  (psi) 


8.50  289.00  0.088 


5.00 

100.00 

0.255 

3. 50 

49.00 

0.520 

3.00 

36.00 

0.708 

2.75 

30.25 

0.  843 

2.25 

20.25 

1.255 

It  should  be  noted  that  the  actual  size  of  the  corner  void  pads  was 
determined  prior  to  calculating  the  static  stress  levels,  due  to  the  physical 
dimension  constraint  on  the  surface  of  the  interior  box.  Should  lower  stress 
levels  be  desired,  the  dimensions  of  the  interior  box  must  be  increased 
accordingly  to  accommodate  the  larger  cushions  which  would  be  required. 
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The  six  temperatures  at  which  drops  were  considered  desirable  were 
-G5,  -20,  20,  70,  110,  and  160°F.  These  are  the  temperatures  specified  by 
the  U.  S.  Army  Natick  Laboratories  for  use  in  testing  bulk  cushioning 
materials.  The  -20°F  temperature  is  especially  important  since  this  is  the 
region  where  some  bulk  cushions  crystallize.  The  four  standard  drop  heights 
utilized  were  12,  18,  24,  and  30  in.  Three  replicates  were  performed  for 
each  set  of  experimental  conditions  with  four  containers  for  each  replicate. 

One  set  of  12  containers  was  utilized  for  temperatures  -20,  70,  and 
160°F,  and  the  other  12  containers  for  20,  110,  and  -85°F.  This  prevents 
any  one  container  system  from  experiencing  the  entire  temperature  range  of 
225°F,  which  is  unlikely  to  occur  in  an  actual  situation.  The  remaining  con¬ 
tainer  system  was  utilized  as  the  prototype  in  the  development  of  the  container 
design. 

Prior  to  conditioning,  the  container  systems  were  instrumi  ted  with 
three  accelerometers  in  the  interior  box,  and  three  in  the  outer  container. 

The  complete  container  systems  were  conditioned  in  environmental  chambers 
(Figure  8)  at  the  required  temperature  for  24  hours  prio^  testing. 

All  tests  were  conducted  at  the  U.  S.  Army  Missile  Research  and 
Development  Command’s  Dynamic  Test  Facility.  Figure  9  shows  the  drop 
tester  and  the  container  positioned  ready  to  drop.  Figure  10  shows  the  test 
apparatus  after  a  drop.  Note  that  the  outer  container  is  almost  as  wide  as  the 
drop  tester  platform.  Therefore,  the  choice  of  size  for  the  outer  container 


Figure  8.  Container  systems  in  environmental  chambers 


Figure  9.  Drop  tester  and  container  system  prior  to  drop, 
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was  constrained.  This,  in  turn,  limited  the  choice  of  cushion  thicknesses  and 

\ 

and  the  ability  to  test  various  thicknesses. 


Experimental  Design 

The  order  of  the  drops  was  randomized  as  much  as  physically  possible. 
One  constraint  on  randomization  was  the  use  of  the  environmental  chambers  for 
temperature  conditioning.  Chamber  space  necessitated  the  drop  test  procedure 
to  consider  all  of  one  temperature  simultaneously.  That  is,  all  three  repli¬ 
cates  for  one  temperature  were  conditioned  as  a  group.  Since  each  box 
required  a  cable  to  be  attached,  just  prior  to  testing  but  after  rem:  /al  from 
the  chamber,  it  was  not  feasible  to  randomize  the  boxes  within  the  replicate. 
However,  the  static  stresses  and  drop  heights  were  randomized  within  each 


replicate.  The  randomized  order  of  the  drops  was  determined  by  the  computer 
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generated  form  shown  in  Figure  11.  One  page  of  this  form  gives  the  order  of 
drops  for  one  replicate  at  a  particular  temperature  on  each  of  four  boxes. 

The  order  of  drops  for  each  box  is  given  in  the  column  under  the  box  number. 
The  letter  for  each  drop  indicates  the  face  (or  stress  level)  on  which  the  box 
is  being  dropped,  while  the  inches  following  the  letter  indicate  the  drop  height. 

The  time  required  to  perform  the  necessary  drops  on  all  six  surfaces 
of  one  box  was  less  than  2  minutes.  Since  temperature  is  involved,  experi¬ 
mentation  was  performed  to  determine  if  a  significant  change  in  temperature 
occurred  prior  to  completion  of  the  sequence  of  six  drops  on  one  container. 
Fortunately,  the  change  in  temperature  was  not  significant  in  the  required 
2-minute  interval. 

The  peak  accelerations  encountered  by  the  interior  box  and  the  total 
box  during  the  432  drops  were  compiled.  The  peak  accelerations  (in  G’s)  for 
the  interior  box  are  given  in  Appendix  A.  Appendix  B  contains  the  peak 
accelerations  for  the  total  box.  It  is  these  basic  test  data  which  will  now  be 
utilized  in  the  development  of  a  mathematical  model  for  the  confined  cushioning 
material  situation.  However,  prior  to  model  development  initiation,  it  is 
necessary  to  consider  the  theoretical  aspects  of  the  situation  which  is  to  be 


modeled, 


MATERIAL:  MINICEL 
TEMPERATURE : 


THICKNESS : 
REPLICATION 
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BOX-1 


C-12" 


F-12" 


E-30" 


B-  2  4  " 


D-30" 


A- 30" 


DROP  HEIGHT 


BOX-2 


BOX- 3 


F-18" 


D-18" 


B-12" 


E-12  " 


C-30" 


B-18" 


A- 1 8  " 


F-30" 


D-24  " 


E- 1 8  " 


A-24"  .*  C-18"  : 


:  BOX- 4  : 


E-24" 


D-12" 


F-24  " 


A-12" 


B-30" 


C-24"  j 


Figure  11.  Randomized  drop  test  sequence 


CHAPTER  IV.  THEORETICAL  CONSIDERATIONS 


A  container  system  is  composed  of  three  elements;  the  exterior  ship¬ 
ping  container,  the  cushioning  material,  and  the  protected  item.  Figure  12A 
depicts  the  system  with  the  three  elements.  By  substituting  springs  for  the 
cushioning  material,  the  system  may  be  represented  as  in  Figure  12B.  To 
study  the  reaction  in  only  one  direction,  all  the  springs  may  be  represented  as 
one  as  in  Figure  12C.  When  a  container  is  dropped,  it  obeys  the  laws  of 
Galileo  and  Newton  by  uniformly  increasing  its  speed  of  descent  under  the  pull 
of  gravity  until  it  strikes  whatever  floor  or  platform  intervenes  at  a  speed 


Figure  12.  Schematics  of  standard  container  system  [3]. 


dependent  upon  the  drop  height  (8. 5  mph  for  a  30-in.  drop  height).  As  shown 
in  Figure  13,  the  container  stops  but  the  protected  item,  also  moving  at 


8.5  mph,  pushes  into  the  cushioning  until  the  resistance  of  the  compressed 
cushioning  is  equal  to  the  force  of  the  moving  item.  The  protected  item, 
slowed  to  a  stop  after  expending  its  energy  of  motion  (derived  from  the  free 
fall)  in  compressing  the  cushioning,  is  returned  to  its  original  position  by  the 
energy  stored  in  the  compressed  cushion. 


START  PACKAGE  PRODUCT  SPRING 

OF  HITS  STOPS  RETURNS 

DROP  FLOOR  SPRING  PRODUCT 

COMPRESSED 


Figure  13.  Free  fall  of  container  system  shown  in  various  stages  [3J. 


The  cushion’s  function  is  to  decelerate  the  protected  item  in  a  few 
inches  which  is  a  fraction  of  the  drop  height  originally  available  for  the  package 
to  accelerate  to  its  impact  velocity.  Since  the  initial  acceleration  was  due  to 
the  force  of  gravity,  the  deceleration  must  be  proportionately  greater,  since 
the  distance  available  is  much  less. 
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The  protected  item,  during  its  deceleration,  experiences  the  same 
pressure  as  the  cushion  which  is  resisting  the  downward  pressure  of  the 
descending  protected  item.  This  force  starts  at  zero  at  the  time  contact  and 
initial  compression  are  made  and  increases  to  a  maximum  when  the  protected 
item  is  stopped  and  returns  to  zero  as  the  compressed  cushioning  material 
returns  the  product  to  its  original  position. 

In  developing  the  impact  response  of  the  protected  item  several 
parameters  were  considered.  From  previous  research,  the  peak  acceleration 
is  related  to  other  parameters  as 

G  =  f(ag,  T,  0,  h) 

where 

G  =  acceleration  (G’s) 

=  static  stress  (psi) 

T  =  thickness  of  cushion  (in.) 

0  =  cushion  temperature  (°F) 

h  =  drop  height  (in.). 

The  objective  of  the  derivation  which  follows  is  to  determine  the  rela¬ 
tionship  of  the  previously  mentioned  parameters  in  die  impact  response  model 
for  the  interior  box.  However,  some  basic  concepts  will  be  discussed  first. 

The  theoretical  model  to  predict  impact  response  is  developed  from  the 
constitutive  relationships  taken  from  viscoelastic  theory.  The  approach  is  to 
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consider  the  viscoelastic  cushioning  material  as  a  standard  linear  solid, 
defined  as  a  material  in  which  the  stress-strain  behavior  can  be  represented 
by  linear  relations  involving  not  only  stress  and  strain  but  also  their  time 
derivatives  of  all  orders.  An  elastic  material  represents  the  special  case  of 
a  linear  material  in  which  the  time  derivatives  can  be  neglected  for  the  range 
of  frequencies  involved.  Materials  which  are  time  dependent  in  their  response 
are  called  viscoelastic. 

In  past  studies  when  a  falling  body  was  dropped  directly  onto  the 
cushion  placed  on  a  rigid  base,  a  Kelvin  viscoelastic  model  was  utilized  to 
depict  this  situation.  The  present  research  is  more  complicated  in  that  the 
falling  body  and  the  cushion  are  packed  together  in  a  wooden  box  which  is 
dropped  onto  the  rigid  base.  The  box  provides  shock  absorbing  or  cushioning 
capacity,  as  does  the  cushioning  material.  The  mechanical  model  used  to 
describe  this  behavior  is  shown  in  F  igure  14.  The  cushioning  material  is 
represented  by  a  Kelvin  model.  The  wooden  box  is  represented  by  a  group  of 
springs  in  parallel  in  which  the  springs  impact  the  rigid  plane.  The  shock  is 


Figure  14.  Interior  box  cushioned  under  impact  loads. 
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transmitted  through  the  springs  representing  the  box  into  the  Kelvin  model 

representing  the  foamed  cushion  material  and  finally  into  the  protected  item. 

Consider  a  container  system  which  is  dropped  a  distance  h .  Inside 

the  container  is  another  box  with  weight  W  and  mass  m.  The  interior  box 

is  cushioned  by  the  cushioning  material  which  is  in  a  tailored  corner  void 

configuration  of  thickness  T .  The  static  stress  between  the  cushion,  with 

total  area  per  side,  A,  and  the  interior  box,  is  represented  by  cr  . 

s 

The  spring  and  dashpot  configuration  on  the  left  in  Figure  14  represents 
a  standard  linear  solid.  The  stiffness  of  the  outer  box  is  represented  by 
parallel  springs  whose  sum  is  k  .  The  relaxation  modulus  according  to 

G 

Flugge  [14]  for  the  standard  linear  solid  maybe  expressed  as 


Vke  +  V 


_t(ke+kl)/7J 


(IV-1) 


where  t  represents  time  and  n  represents  the  viscosity  of  the  cushioning 
material. 

The  equation  oi  motion  of  the  interior  box  after  contact  with  the  cushion 
may  be  expressed  as 


m*x  =  W  -  a  A 
s 


(IV-2) 


where  x  is  the  position  of  the  interior  box  with  respect  to  the  top  of  the 


cushioning  material. 
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If  the  model  used  is  that  of  a  linear  viscoelastic  solid,  then 


a  =  E  (t)€ 
s  rv  ' 


The  strain  e  is  defined  as  the  ratio  of  displacement  x  to  thickness  T  of  the 
cushioning  material  so  that  the  static  stress  may  be  expressed  as 


a 

s 


Er(t)x 


and  the  equation  of  motion  becomes 
AE  (t)x 

mx  =  W - * -  .  (IV-3) 

The  correspondence  principle  described  by  Flugge  [14]  states  that  under 
impact  conditions,  the  modulus  of  elasticity  (which  is  not  time  dependent)  may 
be  substituted  for  the  viscoelastic  relaxation  modulus  in  linear  solids,  giving 

mx  =  W  -  ,  (IV-4) 

where  E  is  the  elastic  counterpart  corresponding  to  Er(t). 

The  general  solution  of  Equation  (IV -4)  may  be  expressed  as 

WT 

x(t)  =  - —  +  A  sin  wt  +  B  cos  wt 
'  '  AE 


elasticity  so  that 


X(t) 


/  2AEr(t)h 


SN  1  + 


WT 


sin  (wt  -  <p) 


Since  only  the  peak  accelerations  are  of  interest  in  this  research,  the 
oscillatory  component  may  be  ignored.  Using  the  series  approximation 


r- -  ,1  1  2  1  3 

’ ^ “  1  +  5X-5X  +I5X  +--- 


the  peak  acceleration  may  be  expressed  as 


*x(t)  =  -g 


hAEr(t)  1(hAEr(t)Y 
1+  WT  ~  2  \  WT  /  +  “  * 


(IV-7) 


Past  research  [10,  17]  has  shown  that  only  the  quadratic  terms  are  signifi¬ 
cant.  The  stiffness  factor  k  in  the  relaxation  modulus  Er(t)  should  be 
many  times  greater  than  the  stiffness  factor  k^  .  Then,  for  impact  loads,  the 
magnitude  of  Er(t)  may  be  approximated  as 


IE  |  «  k  e 
r  e 


C/T] 


1 


‘m 


where  c  is  a  constant.  Equation  (IV-7)  may  then  be  expressed  as 


x(t)  =  -g 


1  + 


hAk  e 
e 


c/77 


,  . hAk  e 
1  I  e 


c/7]N 


WT 


WT 


(IV-8) 


The  quantity  W/A  is  the  static  stress  crg  .  If  the  stiffness  factor  is 


assumed  to  be  approximately  equal  to  the  square  of  the  intensity  of  the  load 


as  measured  by  the  static  stress,  Equation  (IV-8)  may  be  expressed  as 


x(t)  =  -g 


1  + 


ha  e 
s 


c/rj 


,  ,  ha  e 
1  s 


c/V 


2  \  T 


(rv-9) 


For  values  of  a  less  than  approximately  1. 5,  a  transformation  function  of 
s 

the  form  (1  -  cos  a  )  permits  Equation  (IV-9)  to  be  expressed  approximately 
s 

as 


x(t)  *  -  g 


1  + 


he 


c/77 


T  (l-cosffs) 


1  h2 


XU  /.  «  O  2c/77 

-  2  ^2  (!  -cos  *sr  *  e  +  •  •  • 


(IV-10) 


If  a  series  approximation  is  used  for  in  the  form 


x  x2 

e  =  l  +  x  +  —  + 


and  the  viscosity  17  is  inversely  proportional  to  the  temperature  0  ,  then 


x(t)  =  -g 


h(  1  -  cos  a  ) 


1  + 


(C^  C20  +  C302  +  .  .  .  ) 


-  (1  -  cos  ag)2  (C1  +  C.2e  +  c3e2  +  . . .  )2  + 


If  only  arbitrary  constants  are  used,  selected  terms  from  the  following  general 


form  may  be  used  to  approximate  the  peak  acceleration: 
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n  /  \i  m 

G  =  C+  (l-cos(Ts)1  X  D.iJ  .  (IV-11) 


Previous  studies  [10,  17]  on  the  impact  response  of  uneonfined  cushioning 


materials  predicted  an  expression  of  the  form 


„  ,  1/2 

n  C  .h  m 

G  =  C+  2  ‘  ,  £  k  *'  .  ( IV— 12) 

1=1  (<rsT  )1/2  j=l  1J 


Expressions  of  the  form  given  in  Equations  (IV-11)  and  (IV-12)  are 
likely  candidates  for  empirically  curve-fitting  cushioning  system  experimental 
data.  Chapter  V  will  utilize  the  previous  expressions  as  a  framework  from 


which  to  develop  a  general  model  for  the  conditions  under  study 


CHAPTER  V.  MODEL  DEVELOPMENT 


As  stated  previously  in  Chapter  IV,  the  relationship  of  static  stress  to 
weight  of  a  cushioned  item  and  total  cushioned  area  is  given  by  the  equation 

W 

a  =  — 
s  A 

Although  the  19-in.  plywood  cubes  utilized  as  the  dummy  items  within  the 
containers  were  constructed  according  to  government  specifications,  their 
individual  weights  varied  enough  to  cause  the  static  stresses  to  vary  among  the 
24  samples.  Table  I  lists  the  weights  of  the  protected  interior  box,  the 
exterior  box,  and  the  total  box.  The  resulting  variation  in  the  respective 
static  stresses  for  each  of  the  six  surfaces  of  each  cube  are  given  in  Table  II. 

The  data  analysis  procedure  was  determined  prior  to  the  experimenta¬ 
tion.  This  procedure  requires  the  experimental  data  to  be  in  replicates  of 
three  for  each  surface  of  the  experimental  boxes,  as  outlined  in  the  experi¬ 
mental  considerations  chapter.  Consequently,  it  was  necessary  to  devise  a 
method  to  transform  the  experimental  data  for  the  interior  box  to  reflect  a 
common  static  stress  level  for  each  respective  side.  This  was  accomplished 
through  the  use  of  the  following  linear  transformation.  If  the  adjusted  acceler¬ 
ation  is  denoted  by  G  ,  the  uniform  static  stress  by  a  ,  the  measured  peak 

s 
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Surface 


Box  No. 

A 

B 

C 

D 

E 

F 

1 

0.082 

0.500 

0.681 

n 

1.210 

2 

0.088 

0.263 

0.537 

0.731 

1.299 

3 

0.088 

0.261 

0.532 

0.724 

0.862 

1.287 

4 

0.084 

0.250 

0.510 

0.694 

0.826 

1.235 

5 

0.089 

0.265 

0.541 

0.736 

0.876 

1.309 

6 

0.087 

0.258 

0.527 

0.717 

0.853 

1.275 

7 

0.085 

0.252 

0.514 

0.700 

0.833 

1.244 

8 

0.088 

0.262 

0.534 

0.727 

0.866 

1.293 

9 

0.083 

0.248 

0.505 

0.688 

0.818 

1.222 

10 

0.085 

0.254 

0.519 

0.707 

0.841 

1.256 

11 

0.086 

0.257 

0.524 

0.714 

0.849 

1.269 

12 

0.088 

0.261 

0.533 

0.726 

0.864 

1.290 

13 

0.083 

0.246 

0.501 

0.682 

0.812 

1.213 

14 

0.083 

0.248 

0.505 

0.6C3 

0.818 

1.222 

15 

0.088 

0.263 

0.536 

0.729 

0.868 

1.296 

16 

0.080 

0.238 

0.486 

0.661 

0.787 

1.176 

17 

.3.035 

0.254 

0.519 

0.707 

0.841 

1.256 

18 

0.088 

0.264 

0.538 

0.733 

0.872 

1.302 

19 

0.085 

0.254 

0.518 

0.705 

0.839 

1.253 

20 

0.084 

0.249 

0.508 

0.691 

0.822 

1.228 

21 

0.086 

0.255 

0.520 

0.708 

0.843 

1.259 

22 

0.089 

0.265 

0.541 

0.736 

0.876 

1.309 

23 

0.083 

0.246 

0.503 

0.684 

0.814 

1.216 

24 

0.084 

0.251 

0.513 

0.198 

0.831 

1.241 
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acceleration  by  A  ,  and  the  actual  static  stress  by  S  ,  then  the  following 
relationship  may  be  used  to  obtain  the  adjusted  peak  acceleration: 


The  original  data  with  the  associated  static  stresses  are  given  in  Appendix  A. 
The  adjusted  accelerations  are  given  in  Appendix  C. 


Preliminary  Data  Analysis 

With  any  experimental  effort,  a  set  of  observations,  assumed  to  be 
taken  under  the  same  conditions,  may  vary  widely  from  other  observations  or 
be  what  is  known  as  an  outlier.  The  experimenter  must  decide  whether  to 
retain  the  outlier  observation  in  his  computations  or  to  discard  it  as  a  faulty 
measurement.  Even  though  the  data  point  is  discarded  from  computation,  it 
is  still  recorded. 

The  argument  for  exclusion  of  outlying  data  is  that  if  the  data  are  good 
one  simply  loses  some  of  the  relevant  information;  in  this  research  effort,  it 
means  discarding  one  of  the  three  replicated  observations.  If  the  discarded 
observation  is  truly  inaccurate-  then  its  inclusion  would  bias  the  results  and 
the  estimate  of  precision  by  an  unknown,  and  generally  unknowable,  amount. 

Consequently,  it  was  determined  that  a  test  for  outliers  would  be  per¬ 
formed.  The  first  step  in  the  outlier  test  is  to  compute  the  sample  variance, 
while  holding  static  stress  constant,  for  each  set  of  three  replications  of  peak 
accelerations  to  find  which  set  has  the  maximum  sample  variance.  For  the 


set  of  observations  having  the  largest  variance,  each  observation  of  the  set  is 
then  tested  individually  as  a  candidate  for  rejection  as  an  outlier  by  using  the 
statistic 


s 

v 


where 

x^  =  an  individual  observation  in  a  set  of  three  replications 
x  =  the  sample  mean  of  the  three  observations 

s  =  an  independent  external  estimate  of  the  standard  deviation  from 
concurrent  data. 


This  statistic  is  based  upon  an  extension  of  the  extreme  studentized  deviate 
from  the  sample  mean  (Nair  Criterion). 

The  set  of  replications  of  G’s  having  the  maximum  sample  variance 
corresponding  to  a  particular  stress  level  is  eliminated  from  the  calculation 
of  s^.  From  the  remaining  sets  of  replications  of  G's,  s^  is  calculated  with 
the  expression 


til 

where  s.2  is  the  sample  variance  of  the  i  set  of  replications  of  G’s  and  n  is 


the  number  of  stress  levels 
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The  values  of  the  statistic  for  each  observation  in  the  set  of  replicates 
being  tested  are  compared  with  the  appropriate  value  from  a  table  of  per¬ 
centage  points  of  the  extreme  studentized  deviate  from  the  sample  mean,  and 
a  point  is  rejected  as  an  outlier  if  t(calculate(l)  >  t{t>Ue)  . 

If  an  observation  is  rejected  in  the  first  iteration  of  the  outlier  test, 
the  set  of  replicates  to  which  it  belonged  is  no  longer  considered  in  further 
calculations,  but  the  procedure  then  moves  to  the  set  of  replicates  with  the 
next  highest  sample  variance  to  check  for  outliers.  Iteration  is  continued 
until  a  set  of  replications  is  checked  and  no  points  are  rejected.  Those  data 
points  determined  to  be  outliers  are  marked  with  an  asterisk  in  Appendix  B. 
Also,  a  computer  code  for  determining  outliers  as  described  is  contained  in 
Appendix  D. 


Interior  Box  Model  Development 

With  the  outlying  data  points  determined  and  isolated  from  the  compu¬ 
tations,  development  was  initiated  on  the  polynomial  regression  model  which 
would  predict  the  maximum  acceleration  the  interior  box  undergoes  upon 
impact  ( inside  a  container)  for  a  given  temperature  and  drop  height.  Initial 
attempts  to  acquire  a  meaningful  model  with  the  experimental  data  in  its 
original  form  proved  fruitless.  Consequently,  a  search  was  initiated  for  a 
data  transformation  for  the  static  stress  variable.  If  y  is  the  predicted 
acceleration  and  x.  represents  the  static  stress,  the  general  form  of  the 
p-order  polynomial  regression  equation  is: 


They  represent  the  specific  mathematical  model  given  for  each  temperature , 
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Table  III.  Interior  Box  EDCC  Equations 


Temperature 

(°F) 

Drop 

Height 

(in.) 

Box  Equation 

-65 

12 

y  =  34.56  -  40.75(1  -  cos  x)  +  24.83(1  -  cos  x)2 

18 

y  =  43.13  -  45.38(1  -  cos  x)  +  29.40(1  -  cos  x)2 

24 

y  =  32.15  2.97(1  -  cos  x)  +  9.08(1  -  cos  x)2 

30 

y  =  46. 9C  -  58.71(1  -  cos  x)  +  60.42(1  -  cos  x)2 

-20 

12 

y  =  42.26  -  110.62(1  -  cos  x)  +  118.69(1  -  cos  x)2 

18 

y  =  43.43  -  74.75(1  -  cos  x)  +  68.70(1  -  cos  x)2 

24 

y  *  64.63  -  170.13(1  -  cos  x)  +  164.84(1  -  cos  x)*’ 

30 

y  =  62.51  -  123.09(1  -  cos  x)  +  120.88(1  -  cos  x)2 

20 

12 

y  =  42.18  -  109.62(1  -  cos  x)  +  108.26(1  -  cos  x)2 

18 

y  =  41.70  -  55.87(1  -  cos  x>  +  48.88(1  -  cos  x)2 

24 

y  =  43.74  -  49.24(1  -  cos  x)  +  36.23(1  -  cos  x)2 

30 

y  *  57.55  -  130.25(1  -  cos  x)  +  151.13(1  -  cos  x)2 

70 

12 

y  =  32.07  -  58.67(1  -  cos  x)  +  61.64(1  -  cos  x}2 

18 

y  =  38.33  -  43.61(1  -  cos  x)  +  30.49(1  -  cos  x)2 

24 

y  3  36.04  -  9.38(1  -  cos  x)  +  0.09(1  -  cos  x)2 

30 

y  =  60.49  -  145.11(1  -  cos  x)  +  160.52(1  -  cos  x)2 

110 

12 

y  =  27.03  -  49.05(1  -  cos  x)  +  47.24(1  -  cos  x)2 

18 

y  =  28.47  -  13.52(1  -  cos  x)  +  7.67(1  -  cos  x)2 

24 

y  =  46.28  -  73.39(1  -  cos  x)  +  65.84(1  -  cos  x)2 

30 

y  -  36.47  -  45.93(1  -  cos  x)  +  64.80(1  -  cos  x)2 

loJ 

12 

y  =  31.97  -  86.81(1  -  cos  x)  +  96.25(1  -  cos  x)2 

18 

y  =  32.30  -  50.53(1  -  cos  x)  +  53.23(1  -  cos  x)2 

24 

y  *  31.17  -  46.78(1  -  cos  x)  63.77(1  -  cos  x)2 

30 

y  3  45.80  -  54.76(1  -  cos  x)  +  51.37(1  -  cos  x)2 

Figure  15  provides  a  representative  sample  of  an  IDCC  plotted  from  an  equa¬ 
tion  from  Table  III. 

The  IDCC’s  developed  are  applicable  only  to  a  specific  temperature 
and  drop  height  with  2-in.  Minicel  bulk  cushioning.  The  next  step  is  the 


1 
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Figure  15.  Interior  box  IDCC  prediction  at  -65°F  and 
a  30-in.  drop  height. 


development  of  the  general  model  for  2-in.  Minicel  which  has  as  independent 

variables,  drop  height  h,  static  stress  a  ,  thickness  of  cushion  T,  and 

s 

temperature  9  ,  and  as  a  dependent  variable,  G,  the  peak  acceleration. 

The  initial  model.  Equation  (IV-11) ,  proved  deficient  in  representing 
the  nonlinear  characteristics  of  cushion  response  for  the  interior  box.  Con¬ 
sequently,  a  modular  modeling  technique  was  utilized  in  which  each  independ¬ 
ent  variable  was  studied  with  its  effect  upon  the  dependent  variable.  Previous 
research  efforts  [10]  have  shown  that  the  peak  acceleration  is  related  to  an 
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exponential  form  of  temperature  0  .  Similarly,  the  drop  height  h  has  been 
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found  to  enter  the  mathematical  model  as  the  square  root  of  that  value,  while 

thickness  T  is  a  negative  exponential.  Note  that  the  data  transformation 

(l-cos  a  )  used  for  static  stress  in  the  IDCC  equations,  is  also  utilized 
s 

here  for  static  stress. 

A  stepwise  regression  procedure  was  utilized  in  acquiring  the  general 
model  expression.  Draper  and  Smith  [  18 J  offer  this  procedure  as  an  improved 
version  of  the  forward  selection  process  for  variable  selection.  In  the  step¬ 
wise  regression  procedure,  the  variables  already  in  the  equation  are  reeval¬ 
uated  at  each  stage.  A  total  of  45  terms  were  examined  in  the  stepwise 
regression  procedure.  These  terms  were  various  combinations  of  variables 
found  by  previous  research  [9,  10,  11]  to  describe  the  behavior  of  bulk 
cushioning.  Table  IV  lists  the  45  terms.  The  variable  combinations  are 
identified  by  an  x  in  the  appropriate  column.  A  zero  value  in  the  coefficient 
column  indicates  a  variable  combination  which  is  not  in  the  developed  model. 

The  interior  box  general  model  may  be  stated  as: 


G  =  C.+  2  hi/2  2  -rrj-r  £  2  C.  cos  <r  )* 


3  .  2 


i= 0  k 


%  T(l/2+k) 


+  Yj  0n  Z  C  n  (1  "  C0S  °  )m  • 

n=l  m=0  s 


From  this  general  model  it  is  necessary  to  select  the  combination  of 


terms  which  best  describes  the  situation  to  be  modeled.  Each  time  a  variable 


is  added  to  the  interior  box  general  model,  it  is  necessary  to  evaluate  the 


resultant  dynamic  cushioning  curves  to  assure  the  proposed  model  is  provid¬ 
ing  the  hypothesized  U-shaped  curves  of  unconfined  testing  which  do  not 


possess  negative  peak  accelerations  and  the  curves  are  distinct.  Obviously, 
many  of  the  proposed  models  are  similar  in  their  predictive  ability.  However, 
based  upon  first-hand  observations  during  the  experimental  phase,  and 
supported  by  statistical  evaluations  at  each  step  of  the  model  development 
effort,  it  was  determined  that  the  best  interior  box  general  model  to  describe 
the  confined  situation  for  the  cushioning  material  under  consideration  was  a 
constant  and  eight  independent  variables.  The  resulting  nine  term  interior 
box  general  model  is: 


flh1/2  A1/2  e3(l-eos<r)2 

G>  -14. 88  +  7.  ^-0.  58  2-^-0.  54  - - 

2  3 

0(1  -  cos  a  )  0(1-  cos  o  ) 

+  3o.99  — 2  +  66.56  - +  4.10  - 

-  10.22  02  (l  -  cos  a  )  -  0.23  03 

s' 

where 


°F  +  460 


The  confined  general  model  has  a  0. 806  correlation  coefficient.  All 


statistical  tests  have  been  performed  at  an  alpha  level  of  0. 10.  The  model 


should  be  used  only  within  the  ranges  of  the  following  independent  vrriables: 
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h  =  drop  height  from  12  to  30  in. 

og  =  static  stress  range  from  0.088  to  1.255  psi 

9  -  temperature  from  -65  to  160°F. 

The  model  will  predict  accurately  within  the  ranges  given  for  the 
previously  mentioned  variables.  Since  data  were  not  available  outside  the 
ranges  of  the  independent  variables,  the  accuracy  of  the  model  to  predict  in 
those  ranges  is  unknown. 

Figure  16  is  a  typical  example  of  the  predictive  ability  of  the  interior 
box  general  model.  It  is  noted  that  the  curve  is  U-shaped  and  possesses 
characteristics  similar  to  the  unconfined  situation.  Further  comparisons 
between  the  confined  and  unconfined  models  will  be  presented  in  Chapter  VII. 


Figure  16.  Interior  box  general  model  prediction  at  -65°F  and 

a  30 -in.  drop  height. 
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Total  Box  Model  Development 

The  manner  in  which  the  experimental  design  acquired  experimental 
data  permitted  the  acquisition  of  total  box  peak  acceleration  observations  at 
no  additional  testing  cost.  Consequently,  three  accelerometers  were  attached 
on  an  internal  surface  of  the  exterior  box  to  acquire  these  peak  accelerations. 

Since  the  interior  box  varied  in  weight,  the  total  weight  of  the  interior 
and  exterior  box  combined  also  varied.  This  total  weight  variation  was  also 
affected  by  the  moderately  varying  weight  between  each  of  the  individual 
outside  boxes.  A  summary  of  the  individual  box  weights,  both  interior  and 
exterior,  is  given  in  Table  I. 

Since  the  total  box  weights  did  vary,  it  was  hypothesized  that  a  total 
box  model  could  be  developed  that  could  predict  the  peak  accelerations  which 
the  exterior  box  encounters  as  a  function  of  total  box  weight.  Furthermore, 
it  is  intuitive  that  as  box  weight  increases,  the  G-level  will  also  increase. 
What  is  not  known  is  the  effect  of  temperature  upon  the  exterior  box  and  the 
increased  total  box  weight.  It  is  anticipated,  however,  that  the  total  box 
model  will  be  linear  in  nature. 

The  total  box  experimental  data  are  given  in  Appendix  B,  with  truly 
anomalous  readings  identified  with  an  asterisk,  based  upon  an  outlier  pro¬ 
cedure.  The  best  fitting  equation  for  a  particular  drop  height  and  temperature 
situation  was  determined  in  the  same  manner  as  that  utilized  for  the  interior 
box.  The  correlation  coefficients  for  the  individual  box  equations  were  in  the 
0.65  to  0. 80  range.  A  complete  set  of  individual  total  box  equations  is  given 


49 


in  Table  V.  The  computer  program  used  to  obtain  the  individual  total  box 
equations  is  given  in  Appendix  E.  It  is  significant  to  note  that  22  of  the  24 
developed  individual  total  box  equations  possess  a  positive  slope,  verifying 
the  earlier  conjecture  of  increasing  peak  acceleration  as  a  function  of  increas¬ 
ing  weight. 

These  equations  will  predict  total  box  peak  accelerations  for  a  par¬ 
ticular  drop  height,  temperature,  and  total  box  weight.  However,  the  pre¬ 
dictive  range  is  limited  to  the  values  of  the  basic  experimental  data.  Figure 
17  is  a  typical  plot  of  a  total  box  individual  equation. 

To  provide  further  compatibility  between  the  interior  and  exterior 
boxes,  a  general  model  is  desirable  which  depicts  the  peak  acceleration 
anticipated  by  the  total  box.  Consequently,  the  data  utilized  to  develop  the 
individual  total  box  equations  is  combined  to  provide  an  experimental  basis 
for  a  general  total  box  model.  The  total  box  model  development  modifies  the 
procedure  utilized  in  the  development  of  the  interior  box  general  model.  In 
particular,  weight  (W)  is  substituted  for  the  term  (l  -  cos  x). 

The  total  box  model  may  be  stated  as: 


From  this  general  model  it  is  possible  to  identify  the  combination  of 
terms  which  best  describes  the  situation  to  be  modeled.  Each  time  a  variable 
is  added  to  the  total  box  general  model,  it  is  necessary  to  evaluate  the 


Table  V.  Individual  Total  Box  Equations 


Temperature 

(°F) 

Drop  Height 
(in.) 

-65 

12 

18 

24 

30 

-20 

12 

18 

24 

30 

20 

12 

18 

24 

30 

70 

12 

18 

24 

30 

110 

12 

18 

24 

30 

160 

12 

18 

24 

30 

Box  Equation 


-239.44  +  5.78  W 
-281.00  +  7.29  W 
103.33  +  1.24  W 
-122.48  +  4.58  W 


551.45  +11.04  W 
92.67  +  1.06  W 
409.48  +  9.23  W 
-36.99  +  3.52  W 


•  = 

30.66  +  2.21 

W  1 

-  = 

-145.90  +  4.98 

w 

.  = 

94.72  +  1.06 

w 

-  = 

-266.22  +  7.76 

w 

r  = 

-208.96  +  6.67 

w 

r  = 

-440.38  +11.38 

w 

-  = 

-272.13  +  9.19 

w 

102.01  +  2.39 

w 

r  = 

-415.67  +  9.08 

w 

r  = 

26.12  +  2.05 

w 

r  = 

749.16  -  9.66 

w 

-  = 

97.08  +  9.92 

w 

300.21  -  3.30  W 

123.22  +  0.51  W 
237.87  +  6.27  W 
113.78  +  1.31  W 


mm 


O  150 


2  145 


<  135- 

£  130 


55  55.5  56  56.5  57  57.5  58  58.5  59  59.5  60  60.5  61  61.5 

BOX  WEIGHT  (lb) 

Figure  17.  Individual  total  box  prediction  at  -65°F  and 
a  30-in.  drop  height. 


resultant  total  box  line  to  assure  the  proposed  model  is  providing  the  hypothe¬ 
sized  positive  sloped  line.  The  positive  sloped  line  is  based  upon  the  concept 
that  force  is  in  direct  proportion  to  weight.  This  leads  to  the  selection  of  the 
following  total  box  model: 


2  1/2  3  1/? 

G  =  -6.3548  +  0.0723  0  h  W  -  0.0092  0hW  . 


The  confined  total  box  model  has  a  0.772  correlation  coefficient. 

As  can  be  seen,  the  previously  mentioned  model  contains  temperature 
6  ,  drop  height  h ,  and  total  box  weight  W  as  the  parameters  which  vary. 
Consequently,  it  is  possible  to  select  a  particular  temperature,  drop  height, 
and  total  box  weight,  and  predict  the  peak  acceleration  for  those  conditions. 


PEAK  ACCELERATIONS  (G’s) 
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Figure  18  is  an  example  of  the  predictive  ability  of  the  total  box  general 
model. 


Figure  18.  Total  box  general  model  prediction  at  -65‘  F  and 
a  30-in.  drop  height. 


The  immediate  task  then  is  to  v  alidate  the  general  models  which  have 


been  developed  in  this  chapter.  This  will  be  accomplished  in  Chapter  VI, 


CHAPTER  VI.  MODEL  VALIDATION 


In  the  previous  chapter  a  basis  for  a  general  model  was  presented 
which  approximates  the  peak  acceleration  upon  an  item  within  a  cleated  ply¬ 
wood  container  cushioned  with  corner  void  pads  of  2-in.  Minicel.  This  model 
had  the  general  form-. 


G  -  c0+  Z  h</2  Z  -—],/«  Z  «J  Z  CIllr,<l-coS(01 


Li  11  Li  (k+1/2^  iiki' 

f=  0  k=0  ]  j=l  i=0  J 


a  9 

y>  n  v>  ,  .  m 

+  >  0  ),  C  (1  -  cos  a  ) 

u  u  mn'  «' 


n=l  m=0 


Out  of  the  45  terms  generated  by  aie  previous  equation,  eight  terms  were 
retained  using  a  stepwise  regression  procedure.  Before  the  general  model  can 
be  accepted  as  adequately  predicting  impact  response,  it  must  be  validated. 

To  provide  data  for  validation  of  the  general  model,  additional  drop 
tests  were  conducted  in  the  same  manner  as  the  original  drop  tests.  These 
tests  were  made  at  the  same  temperatures  and  static  stress  levels  as  pre¬ 
viously,  but  at  only  one  drop  height,  21-in<  Twenty-one  inches  was  not  one  of 
the  original  drop  heights.  The  21-in.  data  were  not  utilized  in  the  general 
model  development  but  are  unique  data  dropped  at  an  intermediate  drop  height. 
These  data  are  given  in  Appendix  F. 
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Identification  of  Validation  Approach 

The  ultimate  test  of  validity  for  any  model  is  to  assess  the  ability  of 
the  model  to  predict  some  parameter,  in  this  case,  impact  response. 

Natrella  [19]  states  that  "many  statisticians,  when  analyzing  an  experiment 
for  the  purpose  of  testing  a  statistical  hypothesis,  find  that  they  prefer  to 
present  results  in  terms  of  the  appropriate  prediction  limit. "  The  problem, 
of  course,  is  that  every  statistical  test  cannot  be  put  in  the  form  of  a  predic¬ 
tion  limit.  In  general,  tests  that  are  direct  tests  of  the  value  of  a  parameter 
of  the  parent  population  can  be  expressed  in  terms  of  prediction  limits. 

The  choice  then  is  between  a  significance  test  which  provides  a  go/no- 
go  decision,  or  a  prediction  limit  approach  which  provides  much  more  informa¬ 
tion.  The  prediction  limit  procedure  contains  information  similar  to  an  OC 
curve,  and  is  intuitively  more  appealing  than  a  test  of  significance.  If  the 
model  value  is  contained  within  the  prediction  limits,  it  is  possible  to  state 
that  the  model  has  been  validated.  The  width  of  the  prediction  interval  is  a 
good  indication  of  the  firmness  of  the  yes/no  conclusion.  A  great  advantage 
of  the  prediction  limit  approach  is  that  the  width  of  the  limit  is  in  the  same 
units  as  the  parameter.  Consequently,  the  information  is  easy  to  compare 
with  other  information  already  obtained.  Thus,  the  prediction  limit  approach 
will  be  utilized  in  an  attempt  to  validate  both  the  interior  box  model  and  the 


total  box  model 
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Interior  Box  Model  Validation 

The  prediction  limit  method  is  based  upon  the  premise  that  the  general 
model  is  valid  if,  for  selected  values  of  static  stress,  the  general  model 
predicts  values  within  the  prediction  limits.  This  method  generates  predicted 
peak  accelerations  for  both  the  individual  dynamic  cushioning  curve  (IDCC) 
and  the  general  model,  using  the  IDCC  to  establish  prediction  limits. 

Mendenhall  and  Scheaffer  [20]  outline  an  appropriate  equation  for  the 
development  of  prediction  limits  for  a  second-order  equation  which  can  be 
modified  for  cue  interior  box  situation  as  follows: 


=  y  ±  t  s  \ll+  [a 


PLIB  y  ~  a/2 


[a]T  (xTx)  1  [a] 


where  PL  =  interior  box  prediction  limit 
IB 


y=  b0+b1  (1  -cos  crg)  +b2  (1  -cos  ag)‘ 


a  = 


x  = 


(1  -  cos  as) 

4 

(1  -  cos  ag)' 


(1  -  cos  a  \  (l  -  cos  a  \ 

V  V  sij 

(1  -  cos  a  \  / 1  -  cos  a  V 
S2/  \  S2/ 

•  • 

•  • 

•  • 

hi  -  cos  a  \  (l  -  cos  a  \2 

Snj  l  Vj 
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Any  alpha  level  can  be  utilized  in  the  equation  for  the  development  of 
prediction  limits.  However,  to  be  consistent  with  previous  statistical  tests, 
alpha  is  set  at  0.10  for  all  prediction  limit  calculations  in  this  work.  Pre¬ 
diction  limits  are  calculated  at  selected  static  stress  levels  for  the  usual 
combinations  of  temperature  and  a  drop  height  oi  21-in.  The  computer  pro¬ 
gram  for  the  interior  box  validation  is  given  in  Appendix  G. 

Table  VI  provides  a  typical  example  of  the  computer  validation  method 
for  a  temperature  of  20° F.  It  should  be  noted  that  the  static  stress  values  are 
permitted  to  continue  to  1.  538  psi  to  indicate  the  general  form  of  the  curve 
beyond  the  upper  data  limit  of  1. 255.  The  column  entitled  IDCC  represents 
the  predicted  peak  acceleration  based  upon  the  experimental  data  for  a  par¬ 
ticular  temperature  and  a  particular  drop  height.  The  upper  and  lower  pre¬ 
diction  limits  are  developed  around  the  IDCC  function  and  are  shown  in  Columns 
Three  and  Five.  The  general  model  is  then  given  the  specific  temperature  and 
drop  height  under  consideration,  with  the  general  model  output  given  in  Column 
Four.  When  the  general  model  output  is  contained  within  the  prediction  limits, 
the  general  model  is  capable  of  predicting  peak  acceleration  response  with 
sufficient  accuracy  to  be  utilized  in  actual  practice. 

As  can  be  seen  from  Table  VI,  the  values  calculated  by  the  general 
model  fall  within  the  developed  prediction  limits.  In  fact,  at  the  point  of 
maximum  cushioning  (approximately  1.07  psi),  the  IDCC  and  the  model  predict 
within  2.65  G's  of  each  ocher.  Similar  results  are  noted  for  the  other  five 
temperatures  as  shown  in  Table  VII.  In  five  of  the  six  cases,  the  general 
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Table  VI.  Prediction  Limit  Validation  for  Interior  Box 
at  20°  F  and  a  2i-in.  Drop  Height 


Static  Stress 
(psi) 


0.088 


0.288 

0.338 

0.388 

0.438 

0.488 

0.538 

0.588 

0.638 

0.688 

0.738 

0.788 

0.838 

0.888 

0.938 

0.988 

1.038 

1.088 

1.138 

1.188 

1.238 

1.288 

1.338 

1.388 

1.4~8 

1.488 

1.538 


Acceleration  (G) 

IDCC 

Lower -P 

Model 

Upper-P 

49.60 

.  36.03 

45.17 

63.16 

35.52 

jpppfej 

62.48 

48.15 

34.78 

I 

61.53 

47.07 

33.82 

I 

43.40 

60.32 

45.77 

32.65 

42.49 

53.89 

44.28 

31.28 

41.44 

57.28 

42.62 

29.73 

40.27 

55.52 

40.83 

28.01 

39.01 

53.65 

38.94 

26.15 

37.67 

51.72 

36.98 

24.19 

36.27 

49.77 

35.00 

22.17 

34.85 

47.83 

33.03 

20.13 

33.44 

45.93 

31.12 

18.12 

32.05 

44.12 

29.31 

16.20 

30.73 

42.41 

27.65 

14.44 

29.49 

40.85 

26.17 

12.88 

28.38 

39.46 

:-.h.93 

11.59 

27.41 

38.27 

23.96 

10.61 

26.62 

37.32 

23.31 

9.97 

26.04 

36.66 

23.02 

9.70 

25.68 

36.33 

23.11 

9.80 

25.58 

36.42 

23.63 

10.27 

23.76 

37.00 

24.61 

11.05 

26.24 

38.17 

26.07 

12  .09 

27.04 

40.05 

28.04 

! 

13.33 

28.17 

42.74 

30.52 

14.69 

29.64 

46.35 

33.55 

16.15 

31.47 

50.95 

37.12 

17.67 

33.66 

56.56 

41.24 

19.27 

36.21 

63.20 

j  45.90 

20.95 

39.12 

70.85 

SSs 


Table  '.’ll.  Summary  of  Interior  Box  Prediction  Limit 
Validation  Results 


Temperature 

(°F) 

General  Model  Within 
Prediction  Limits 

Prediction  A  at  Maximum 
Cushioning  Level 
(G's) 

-65 

Yes 

-0.05 

-20 

Yes 

5.34 

20 

Yes 

2.56 

70 

Yes 

6.82 

110 

Yes 

4.12 

160 

Yes 

1.37 

model  predicted  higher  than  the  IDCC,  providing  for  a  conservative  prediction 
of  the  cushioning  ability  of  the  particular  material  under  the  specified  condi¬ 
tions.  The  remaining  case  (-..  °F)  indicated  an  optimistic  prediction  for  the 
model.  In  cushion  design  it  is  desirable  to  predict  in  the  conservative  direc¬ 
tion  rather  than  in  the  optimistic  direction.  In  general,  the  results  presented 
in  Table  VII  are  very  acceptable. 


Total  Box  Model  Validation 

The  u>tal  box  model  validation  is  somewhat  less  complicated  than  the 
interior  box  model  validation.  The  primary  difference  between  the  two  models 
is  that  the  interior  box  model  is  quadratic  while  the  total  box  model  is  linear. 
Consequently,  the  total  box  model  prediction  limits  are  acquired  from  the 


first-order  equation  modified  as  follows: 
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where  PL  =  total  box  prediction  limit 


y  -  Vbiw 


The  alpha  level  utilized  for  total  box  validation  purposes  is  0. 10.  The  com¬ 
puter  program  for  total  box  validation  is  given  in  Appendix  H. 

Table  VIII  provides  a  typical  example  of  the  prediction  limits  for  a 
temperature  of  20°F.  The  weight  range  for  the  total  box  is  from  55. 13  to 
61. 31  lb.  Consequently,  a  weight  range  of  55  to  61. 5  lb,  by  0.  5-lb  increments, 
was  selected  for  the  development  of  prediction  limits. 

As  can  be  seen  in  Table  VIII,  the  individual  total  box  line  (ITBL) 
ranges  from  144  G’s  at  55  lb  to  164  G's  at  61.5  lb.  Over  the  same  weight 
range,  the  total  box  model  predicts  values  of  158  to  177  G’s.  The  slope  of 
both  linear  equations  is  noted  to  be  positive  which  appears  to  be  intuitively 
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Table  VIII.  Prediction  Limit  Validation  for  Total  Box  at  20°F  and 

a  21 -in.  Drop  Height 


Weight 

(lb) 

Acceleration  (G) 

ITBL 

Lower 

Model 

Upper 

55.00 

143.95 

29.68 

157.99 

258.23 

55.50 

145.44 

34.08 

159.48 

256.80 

56.00 

146.93 

38.16 

160.98 

255.71 

56.50 

148.42 

41.89 

162.47 

254.96 

57.00 

149.92 

45.24 

163.96 

254.59 

57.50 

151.41 

48.20 

165.46 

254.62 

58.00 

152.90 

50.74 

166.95 

255.05 

58.50 

154.39 

52.86 

168.45 

255.91 

59.00 

155.88 

54.55 

169.94 

257.21 

59.50 

157.37 

55.80 

171.44 

258.93 

60.00 

158.86 

56.82 

172.93 

261.09 

60.50 

160.35 

57.02 

174.42 

263.68 

61.00 

161.84 

57.01 

175.92 

266.67 

61.50 

163.33 

56.60 

177.41 

270.06 

correct,  since  a  heavier  falling  item  is  expected  to  impact  a  rigid  surface  with 
more  G's  than  a  lighter  falling  item.  It  is  also  noted  that  the  model  values 
fall  within  the  prediction  limits  associated  with  the  ITBL  for  this  particular 
case.  In  a  similar  fashion,  the  model  values  for  the  remaining  five  cases  all 
fall  within  their  respective  prediction  limits. 

Table  IX  further  substantiates  the  intuitive  feeling  expressed  in  the 
previous  paragraph.  The  slope  for  all  six  ITBL  equations  is  positive,  indicat¬ 
ing  that  G’s  increase  as  a  function  of  increased  weight.  Table  X  presents  the 
peak  acceleration  ranges  for  the  ITBL  and  the  model  for  the  six  different 


Table  IX.  Regression  Coefficients  for  ITBL  for 
21-in.  Drop  Height 


Temperature 

(°F) 


Coefficient 

bo 

n 

-195.04 

6.14 

-49.02 

3.30 

-20.00 

2.98 

-235.70 

6.27 

65.21 

1.63 

-51.44 

3.46 

Table  X.  Model  and  ITBL  G-level  Ranges  for 
21-in.  Drop  Height 


Temperature 

(°F) 

ITBL 

Model 

-65 

142.85  -  185.85 

135.60  -  153.66 

-20 

132.43  -  155.52 

149.64  -  169.50 

20 

143.96  -  164.83 

157.99  -  178.91 

70 

108.92  -  152.78 

161.55  -  182.92 

110 

154.67  -  166.05 

157.82  -  178.71 

160 

139.09  -  163.33 

143.46  -  162.52 

temperature  levels.  All  six  cases  are  noted  to  exhibit  the  expected  near 
parallel  line  phenomenon,  that  is,  the  two  lines  are  nearly  parallel  and  do  not 
intersect  over  the  data  range.  Also,  in  five  of  the  six  cases,  the  model  predicts 
higher  peak  accelerations  than  the  ITBL,  resulting  in  conservative  designs. 

The  only  case  in  which  this  is  not  true  is  for  -65°F.  Since  this  case  is  at  one 
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of  the  temperature  extremes,  it  is  possible  that  an  unidentified  physical 
phenomenon  is  taking  place.  However,  a  more  plausible  explanation  is  that 
erratic  observations  due  to  uneven  drops  have  caused  a  distortion  in  the  data 
for  this  one  temperature.  It  should  be  noted  that  a  similar  situation  exists 
for  the  interior  box  at  -65°F,  resulting  in  optimistic  designs. 

Thus,  the  validation  for  the  total  box  model  has  shown  that  all  model 
predictions  are  within  the  specified  prediction  limits.  Hence,  it  is  believed 
that  the  total  box  model  is  an  adequate  means  to  predict  peak  acceleration  for 
intermediate  temperature  and  drop  height  requirements  for  a  total  box 
possessing  similar  weight  and  size  characteristics. 


CHAPTER  VII.  RESEARCH  FINPiNGS 


In  Chapter  V  two  general  iTnpact  response  models  were  developed,  one 
for  the  interior  box  and  one  for  the  total  box.  Validation  oi  these  two  general 
models  was  accomplished  in  Chapter  VI.  This  Chapter  will  demonstrate  the 
areas  in  which  these  two  models  may  be  utilized. 

Integration  of  the  Two  General  Models 
As  previously  noted,  the  interior  b^>x  g~:.eral  model  prediction  of  peak 
acceleration  is  a  function  of  static  stress,  temperature,  drop  height,  and 
cushion  thickness.  The  total  box  general  model  prediction  of  peak  acceleration 
is  a  function  of  weight,  temperature,  and  drop  height.  It  is  desirable  to  view 
the  two  general  models  as  an  integrated  pair,  in  which  the  interior  box  general 
model  is  a  subset  of  the  total  box  general  model.  Then,  it  is  possible  to 
isolate  the  shock  absorption  capabilities  of  the  cushioning  material  by  taking 
the  difference  between  the  total  box  model  results  and  the  interior  box  model 
results  for  selected  values  of  static  stress. 

Since  the  two  models  were  developed  with  different  (but  related) 
parameters,  to  achieve  compatibility  it  is  necessary  to  utilize  the  basic 


equation: 


64 


W 

°s  =  A 

where 

a  =  an  interior  box  static  stress  value 
s 

W  =  weight  of  the  total  box 
A  =  footprint  of  the  tot.u  box. 

The  footprint  of  the  total  box  is  the  surface  area  of  that  portion  which  makes 
contact  with  the  rigid  surface.  The  experimental  total  boxes  possessed  sur¬ 
face  areas  of  172. 5  in.2.  Since  the  weight  range  of  the  experimental  boxes 
was  from  55  to  61. 5  lb,  substitution  of  these  values  in  the  previous  equation 
gives  the  range  of  static  stress  values  for  comparison  purposes. 

Utilizing  the  21-in.  validation  data  as  an  example  (Tables  XI  through 
XVI) ,  it  is  observed  that  the  interior  box  has  not  reached  the  optimum  loading 
point  at  any  temperature  level.  The  optimum  point  would  be  an  inflection  point 
on  the  curve  plotted  from  the  interior  box  column.  In  fact,  the  interior  box 
peak  accelerations  are  decreasing  slightly  with  increases  in  total  box  weight. 
Furthermore,  in  each  case  the  cushion  still  possesses  sufficient  ability  to 
absorb  most  of  the  increase  in  peak  acceleration  which  occurs  as  a  function  of 
increased  weight.  As  expected,  the  total  box  peak  accelerations  continue  to 
increase  as  a  function  of  increased  total  box  weight.  The  computer  code 
utilized  to  develop  Tables  XI  through  XVI  is  contained  in  Appendix  I. 


Table  XI.  Integrated  Confined  Model  Data  for  -65°F  and 
a  21-in.  Drop  Height 


Box  Weight 
(it) 

Acceleration  (G)  j 

Total  Box 

Interior  Box 

Cushion 

■a 

135.59 

41.38 

94.20 

Ipspi 

136.88 

41.33 

95.55 

HH 

138.17 

41.27 

96.90 

56.50 

139.46 

41.21 

98.25 

57.00 

140.75 

41.15 

99.60 

57.50 

142.04 

41.09 

100.94 

58.00 

143.33 

41.03 

102.29 

58.50 

144.62 

40.97 

103.64 

59.00 

145.91 

40.91 

105.00 

59.50 

147.20 

40.85 

106.35 

60.00 

148.49 

40.79 

107.7  0 

60.50 

149.78 

40.73 

109.05 

61.00 

151.07 

40.67 

110.40 

61.50 

152.37 

40.61 

111.75 

Table  XII.  Integrated  Confined  Model  Data  for  -20° F  and 
a  21-in.  Drop  Height 


Box  Weight 
(lb) 

Acceleration  (G)  f 

Total  Box 

Interior  Box 

Cushion 

■i 

149.64 

42.27 

107.37 

mmm 

151.06 

42.20 

108.85 

WSSM 

152.47 

42.14 

110.33 

56.50 

153.89 

42.08 

111.81 

57.00 

155.31 

42.02 

113.29 

57.50 

156.73 

41.96 

114.77 

58.00 

158.15 

41.89 

116.25 

58.50 

159.57 

41.83 

117.73 

59.00 

160.98 

41.77 

119.21 

59.50 

162.40 

40.70 

120.70 

60.00 

163.82 

41.64 

122.18 

60.50 

165.24 

41.57 

123.66 

61.00 

166.66 

41.51 

125.15 

61.50 

168.07 

41.44 

126.63 

in?  '  ™rf<gS&5f?BS 


Table  XIII.  Integrated  Confined  Model  Data  for  20°  F  and 
a  21-in.  Drop  Height 


Box  Weight 
(lb) 


i5.00 

>5.50 

16.00 

56.50 
57.00 

57.50 
58.  On 
>8.50 
>9.00 
>9.50 
>0.00 
>0.50 
>1.00 
>1.50 


Acceleration  (G) 

Total  Box  Interior  Box  Cushion 

157.99  41.85  116.13 

159.48  41.79  117.69 

160.98  41.73  119.24 

162.47  41.  t/  120.80 

163.96  41.60  122.36 

165.46  41.54  123.92 

166.95  .  41.47  125.47 

168.45  41.41  127.03 

169.94  41.35  128.59 

171.44  41.28  130.15 

172.93  41.21  131.71 

174.42  41.15  133.27 

175.92  41.08  134.83 

177.41  41.02  136.39 


Table  XIV.  Integrated  Confined  Model  Data  for  70°F  and 
a  21-in.  Drop  Height 


Box  Weight 
(lb) 

55.00 

55.50 
56.00 

56.50 
57.00 

57.50 
58.00 

58.50 
59.00 

59.50 
60.00 

60.50 
61.00 

61.50 


|  Acceleration  (G)  ! 

Total  Box 

Interior  Box 

Cushion 

161.55 

39.65 

121.89 

163.07 

39.59 

123.48 

164.60 

39.54 

125.06 

166.13 

39.48 

126.65 

167.65 

39.42 

128.23 

169.18 

39.36 

129.82 

170.71 

39.30 

131.  U 

172.23 

39.23 

132.99 

173.76 

39.17 

134.58 

175.29 

39.11 

136.17 

176.81 

39.05 

137.76 

178.34 

38.99 

139.35 

179.86 

38.92 

140.94 

181.39 

38.86 

142.52 
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Table  XV.  Integrated  Confined  Model  Data  for  110°  F  and 
a  21 -in.  Drop  Height 


Box  Weight 
(lb) 

Acceleration  (G) 

Total  Box 

Interior  Box 

Cushion 

55.00 

157.81 

36.47 

121.33 

55.50 

159.30 

36.42 

122.88 

56.00 

160.80 

36.37 

124.42 

56.50 

162.29 

36.32 

125.97 

57.00 

163.78 

36.26 

127.51 

57.50 

165.27 

36.21 

129.06 

58.00 

166.77 

36.16 

130.60 

58.50 

168.26 

36.10 

132.15 

59.00 

169.75 

36.05 

133.70 

59.50 

171.24 

36.00 

135.24 

60.00 

172.74 

35.94 

136.79 

60.50 

174.23 

35.88 

138.34 

61.00 

175.72 

35.83 

139.89 

61.50 

177.21 

35.77 

141.43 

Table  XVI.  Integrated  Confined  Model  Data  for  160°F  and 
a  21-in.  Drop  Height 


Box  Weight 
(lb) 

Acceleration  (G) 

Total  Box 

Interior  Box 

Cushion 

55.00 

143.45 

112.82 

55.50 

144.81 

114.22 

56.00 

146.18 

30.55 

115.62 

56.50 

147.54 

30.51 

117.02 

57.00 

148.90 

30.48 

118.42 

57.50 

150.26 

30.44 

119.82 

58.00 

151.62 

30.40 

121.22 

58.50 

152.99 

30.36 

112.62 

59.00 

154.35 

30.32 

124.02 

59.50 

155.71 

30.28 

125.42 

60.00 

157.07 

30.24 

126.83 

60.50 

158.43 

30.20 

128.23 

61.00 

159.79 

30.16 

129.63 

61.50 

161.16 

30.12 

131.03 

68 


Temperature  Effects 

Perhaps  the  most  significant  finding  of  this  research  effort  is  the 
column  entitled  "Cushion”  in  Tables  XI  through  XVI.  Heretofore,  the  actual 
peak  acceleration  absorbed  by  a  cushion  as  a  function  of  weight  and  tempera¬ 
ture  has  never  been  calculated.  Consequently,  Figure  19  illustrates  the 
effect  of  temperature  upon  the  interior  box  cushion  as  a  function  of  total  box 
weight.  It  should  be  noted  that  the  Minicel  cushion  performs  well  at  the 
lowest  temperature,  -65°F.  This  low  temperature  phenomenon  is  not  common 
to  all  cushioning  materials,  but  indicates  some  unique  characteristics  for  the 
Minicel  material.  Further  observation  indicates  that  the  Minicel  material 
performs  better,  from  a  cushioning  standpoint,  at  approximately  70° F. 

This  phenomenon  may  be  due  to  the  closed  cell  construction  of  the  foam 
itself.  Minicel  is  composed  of  tiny  closed  cells  in  which  air  is  entrapped. 
Compressing  the  closed  cell  is  comparable  to  compressing  a  balloon.  Two 
fact/'  j  in  the  construction  of  the  cell  which  contribute  to  the  cushioning 
ability  of  the  foam  material  are  the  entrapped  air  and  the  walls  of  the  cell. 

As  the  foam  material  is  cooled,  the  air  contracts  and  causes  the  cell  to  be 
compressed  even  though  the  cell  walls  are  somewhat  rigid  from  the  cold 
temperature.  At  the  higher  temperatures,  the  walls  of  the  cell  become  softer 
and  are  able  to  be  compressed  with  less  force  than  at  lower  temperatures. 

At  the  very  low  temperatures  the  contraction  of  air  is  a  factor  and  at  higher 
temperatures  the  flexibility  of  the  cellular  walls  is  a  cushioning  factor  [21], 


PEAK  ACCELERATION  (G’s)  ABSORBED 
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Figure  20  illustrates  the  effect  of  temperature  upon  the  interior  box 
peak  accelerations.  It  is  evident  that  the  total  box  weight  has  a  much  smaller 
effect  than  temperature.  Furthermore,  the  interior  box  experiences  the 
greatest  peak  acceleration  near  -20°F,  and  the  smallest  peak  accelerations 
at  the  high  temperature  extreme,  160°F. 

Figure  21  depicts  the  total  box  peak  accelerations  as  a  function  of 
temperature.  In  this  situation,  the  outside  container  is  seen  to  be  affected  by 
the  different  temperature  levels,  which  means  the  wooden  construction  serves 
as  a  better  shock  absorber  at  the  cold  and  hot  extremes  than  at  ambient 
temperature.  The  wood  is  considered  to  be  a  closed  cell  composite  [20  J . 

As  a  closed  cell  material  it  undergoes  changes  as  the  temperature  varies. 

One  additional  effect  of  temperature  is  observed  when  the  data  for  a 
particular  item  weight  from  Tables  XI  through  XVI  is  selected.  A  summary 
of  these  data  is  provided  in  Table  XVII  for  a  55-lb  item.  The  development  of 
an  additional  column  in  Table  XVII  which  gives  the  percent  of  total  box  GTs 
absorbed  by  the  cushion  reveals  that  the  Minicel  cushion  performs  best  at  the 
highest  temperature,  160°F.  It  is  noted  that  at  the  low  temperature,  -65°F, 
the  cushion  absorbs  69%  of  the  total  box  G’s  available  increasing  slowly  until 
a  level  of  79%  is  achieved  at  160° F.  Hence,  the  cushion  absorbs  considerably 
more  of  the  available  shock  at  the  higher  temperature  than  at  the  lower  tem¬ 
peratures.  Similar  results  in  cushion  absorption  occur  for  varying  item 


weights 


TEMPERATURE  (°F) 


Figure  20.  Interior  box  peak  accelerations  as  a  function  of 
temperature  and  trtal  box  weight  for  a  21-in.  drop  height. 


ERATION  (G's) 
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Table  XVII.  Cushion  Absorption  Percentage  as  a  Function  of  Temperature 
for  a  Selected  Item  Weight  and  a  21-in.  Drop  Height 


Temperature 

Acceleration 

(G) 

Percent  of 
G's  Absorbed 

Total 

Box 

Interior 

Box 

Absorbed  by 
Cushion 

-65 

136 

41 

94 

69 

-20 

150 

42 

107 

71 

20 

158 

42 

116 

73 

70 

162 

40 

122 

75 

110 

158 

36 

121 

77 

160 

. 

143 

31 

113 

79 

Confined  Versus  Unconfined  Comparisons 
It  is  generally  accepted  that  designing  cushioning  systems  from  uncon¬ 
fined  (flat  pad)  drop  tests  will  result  in  conservative  (too  much  cushion) 
designs.  The  problem  has  been  that  the  magnitude  of  conservatism  which  was 
experienced  was  unknown.  Consequently,  knowledgeable  cushioning  system 
designers  continue  to  utilize  the  best  source  currently  available  (i.e.,  uncon¬ 
fined  data) . 

The  results  of  this  research  permit  a  comparison  between  confined  and 
unconfined  test  results.  As  previously  noted.  McDaniel  [10]  developed  a 
general  Minicel  cushioning  model  based  upon  unconfined  data.  Both  this 
research  and  McDaniel's  utilized  Minicel  cushioning  material,  with  each 


including  the  2-in.  thickness.  Since  both  models  have  generalized  the 


temperature  and  drop  height  parameters,  it  is  possible  to  compare  the  uncon¬ 
fined  interior  box  general  model  results  with  McDaniel's  unconfined  general 
model  results. 

Typical  results  are  shown  in  Figures  22  through  24,  where  it  is  noted 
that  the  unconfined  model  predicts  peak  accelerations  above  the  confined  model 


Figure  22.  Comparison  of  confined  and  unconfined  peak 
accelerations  at  -65°F  with  a  30-in.  drop  height. 
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Figure  23.  Comparison  of  confined  and  unconfined  peak 
accelerations  at  70°F  with  a  30-in.  drop  height. 


when  one  considers  the  optimum  cushioning  point  as  the  confined  model  mini¬ 
mum.  This  confirms  the  original  hypothesis  that  unconfined  drop  tests  result 


in  conservative  designs.  It  is  further  noted  in  Figures  22  through  24  that  the 
optimum  cushioning  point  (lowest  peak  acceleration)  does  not  occur  at  the 
same  static  stress  level  for  the  unconfined  and  confined  cases.  In  fact,  when 
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Figure  24.  Comparison  of  confined  and  unconfined  peak 
accelerations  at  160°F  with  a  30-in.  drop  height. 

the  optimum  cushioning  point  for  the  confined  case  is  to  the  left  of  the  uncon¬ 
fined  point  (Figure  22) ,  the  two  curves  intersect  to  the  right  of  both  optimum 
points.  When  the  confined  optimum  cushioning  point  is  to  the  right  of  the 
unconfined  optimum  point,  the  two  curves  intersect  prior  to  the  confined 
optimum  point  (F  igure  23) . 

Table  XVIII  summarizes  the  optimum  cushioning  point  at  the  six 
standard  temperatures  for  the  unconfined  and  confined  cases  for  a  30 -in. 


Table  XVIII.  Summary  of  Uneonfined  Versus  Confined  Peak 
Acceleration  Minimums  for  a  30-in.  Drop  Height 


Unconfined 


Confined 


Temperature 

(°F) 

Static  Stress 

G's 

Static  Stress 

G's 

-65 

1.238 

40 

1.038 

33 

-20 

0.938 

43 

0.988 

33 

20 

0.688 

42 

1.088 

32 

70 

0.538 

40 

1.088 

31 

110 

0.388 

'  39 

1.038 

31 

160 

0.338 

44 

0.938 

30 

rJ 


drt  p  height.  For  five  of  the  six  temperatures,  the  peak  acceleration  minimum 
for  the  confined  case  is  at  a  larger  static  stress  value  than  for  the  unconfined 
case.  Only  at  -65°F  is  the  peak  acceleration  minimum  at  a  lower  static  stress 
value  for  the  confined  case.  Furthermore,  irrespective  of  the  temperature  or 
static  stress  value,  the  peak  acceleration  minimums  for  the  confined  case  are 
always  lower  than  those  for  the  unconfined  case.  Hence,  once  again  the  con¬ 
servative  nature  of  the  unconfined  approach  is  identified. 

One  additional  point  concerning  the  confined  and  unconfined  general 
models  warrants  mentioning.  The  general  models  available  permit  the  cushion 
designer  to  acquire  cushion  design  information  at  any  desired  intermediate 
drop  height  value  between  12  and  30  in.  Similar  selections  may  be  exercised 
for  the  range  of  temperature  between  -65  and  160° F. 


Figure  25  is  an  illustration  of  the  confined  and  unconfined  general 


models’  ability  to  predict  at  the  drop  height  value  of  20-in.  for  a  temperature 
of  110°F.  The  usual  phenomena  between  the  two  general  models  are  noted. 


Figure  25.  Comparison  of  unconfined  and  confined  peak 
accelerations  at  110° F  with  a  21-in.  drop  height. 


CHAPTER  VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  scientific  evolvement  of  cushioning  design  is  beginning  to  reach 
fruition.  Mindlin’s  [1]  original  scientific  basis  continues  to  serve  as  a  frame¬ 
work  for  advances  into  understanding  the  phenomena  of  cushioning.  The  equa¬ 
tions  of  motion  first  derived  by  Mindlin  served  as  a  basis  for  the  work  con¬ 
tained  herein. 

This  research  has  resulted  in  the  development  of  a  general  mathemat¬ 
ical  model  for  a  confined  cushioning  system  and  a  general  mathematical  model 
for  the  exterior  container  which  surrounds  the  confined  corner  void  configured 
cushions. 


Conclusions 

The  objective  of  this  research  was  to  develop  the  methodology  for 
modeling  the  impact  response  for  the  Minicel  cushioning  material  in  the 
unconfined  state.  This  objective  was  satisfied  through  the  development  of  an 
experimental  drop  test  design,  conducting  an  extensive  drop  test  program,  and 
then  modeling  the  resultant  test  data. 

The  developed  general  model  portraying  the  peak  accelerations  of  the 


confined  cushioning  system  utilizing  2  lb/ft3  Minicel  material  is 
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°  -  V  2*,/2  Z  -n^jZ^Z  ‘Wi-coscg1 

u  i=0  k=0  T1  1/  }j=0  i=0  ljic£  s 


+  Yj  d  Lc  (1_  cos  a ) 

“  mn'  s' 

n=l  m=0 


The  model  is  predicated  upon  viscoelastic  theory  and  incorporates  the  effect  of 
drop  height,  static  stress,  thickness,  and  temperature  of  the  cushion  upon  the 
peak  acceleration  of  a  confined  cushioning  system.  This  general  model  may 
provide  the  basic  underlying  structure  of  peak  acceleration  for  any  of  the  other 
bulk  cushioning  materials  utilized  in  the  confined  configuration  for  shock 
isolation.  Utilizing  the  developed  confined  Minicel  model,  it  is  possible  to 
nredict  the  peak  accelerations  which  the  protected  item  will  experience.  As 
expected,  considerable  disparity  exists  between  the  confined  and  unconfined 
(McDaniel’s  research)  results  for  a  protected  item.  It  has  been  known  for 
some  time  that  the  unconfined  data  result  in  conservative  cushioning  designs. 
These  confined  results  provide  a  measure  of  this  conservatism. 

In  addition,  a  general  model  was  developed  portraying  the  peak  accel¬ 
erations  experienced  by  the  total  container.  This  model  may  be  stated  as 


v  i  hf/2£  •*  i  <v<w)* 

i=  0  j=l  i=0  J 


This  model  incorporates  the  effects  of  drop  height,  temperature,  and  total  box 


weight  upon  the  peak  acceleration  of  the  lulal  container. 


The  total  container  model  places  the  temperature  concept  in  perspec¬ 
tive,  since  it  is  possible  to  identify  the  cushion  absorption  capability  as  a 
function  of  changing  temperatures.  Interestingly,  the  cushion  is  most  effec¬ 
tive  at  the  high  temperature  extreme,  which  is  a  unique  characteristic  of 
polyethylene  cushions. 

It  is  believed  that  the  combination  of  the  two  developed  models  provides 
the  cushioning  system  designer  with  another  measure  of  cushion  performance 
to  assist  in  the  design  of  cushioning  systems.  In  particular,  the  assessment 
of  temperature  controlled  confined  and  unconfined  cushions  provides  informa¬ 
tion  unavailable  until  this  time. 

Recommendations 

Based  upon  the  consistent  results  obtained  in  this  investigation,  it  is 
recommended  that  additional  drop  tests  be  performed  to  extend  these  results 
to  other  bulk  cushioning  materials.  This  would  provide  the  cushioning  system 
designer  with  a  more  valid  indication  of  the  materials*  performance  in  the 
confined  state. 

It  also  appears  reasonable  to  develop  an  expanded  experimental  design 
which  includes  containers  from  the  10-lb  range  up  to  200  lb.  This  would 
permit  the  extension  of  the  total,  box  concept  to  the  entire  static  stress  range 
ordinarily  desired.  Perhaps  the  experimental  design  could  emphasize  data 
compaction  between  varying  box  sizes  as  a  means  to  reduce  the  voluminous 


data  problem. 
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The  testing  of  containers  in  a  controlled  temperature  environment  is 
very  costly.  However,  based  upon  the  significant  temperature  results  obtained 
by  McDaniel,  and  the  drop  tests  conducted  in  this  research,  it  is  recommended 
that  any  additional  confined  cushioning  testing  include  the  minimum  six  tem¬ 
peratures  utilized  in  this  work. 

It  is  also  possible  that  the  general  methodology  provided  in  this 
research  effort  may  be  applicable  to  a  relatively  recent  approach  to  cushion¬ 
ing,  which  involves  a  material  in  the  form  of  bubbles  in  which  air  is  entrapped. 
Since  little  is  known  about  these  bubbled  cushioning  materials'  performance  in 
the  unconfined  state,  perhaps  it  would  be  best  to  progress  to  the  confined  state 
immediately.  The  methodology  is  now  available  to  make  this  analysis,  if  the 
cushioning  system  designer  desires. 
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APPENDIX  A.  INTERIOR  BOX  PEAK  ACCELERATIONS 
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Table  A-5.  Original  Peak  Accelerations  with  Individual  Static  Stresses  for  Interior  Box  at  110°F 
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APPENDIX  B.  TOTAL  BOX  PEAK  ACCELERATIONS 
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Table  B-6.  Individual  Weights  and  Peak  Accelerations  of  Total  Boxes  at  160°F 
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Table  C-2.  Peak  Accelerations  After  Conversion  to  Common  Static  Stress  for  Interior  Box  at  -20°F 
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Table  C-6.  Peak  Accelerations  After  Conversion  to  Common  Static  Stress  for  Interior  Box  at  160°F 
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,.«.*♦*********+*♦**************************♦** ***************  *********** 

E  InTeRiOR  oOX  DATA  ANALYSIS 

e 

t ********** ************************** **********************  ************ 

e 

e 

C  1'iAili  DRIVER  ReaUS  THE  CUSHION  MATERIAL  DATA  RECORD  CO  LISTING  OF 

e  temperature*  drop  height*  Thickness*  stress-level*  g-jalue* 

e  KePLaCAT ION  MNP  MATtRI«L  TYPE  (In  THIS  ORDER). 

E  INITIALIZES  ARRaYS  CALl  OuTlR  (OUTLIER  SUBPROGRAM)*  A  <D  CALLS 

L  CVkE6  (CURVILINEAR  REGRESSION  SUbPHOGRAM). 

c 

c  iHt  output  is  in  the  form  of: 

l  for  a  particular 

L  1.  DROP  HEIGHT* TEMPERATURE* AND  MATERIAL  TttlCKNESa  A 

C  TAoLE  eUNTAINInG  The  STRESS  LEVELS  AND  G-VALUiS* 

e  And  AnY  POINT  THAT  »AS  REJECTED  BY  OUTuR  IS  A-SO  LISTED* 

E 

E  2.  The  F  -  STATISTIC  anD  THe  FIRST*  SECOND*  THIRu  AND 

o  FOURTH  uEGREE  POLYNOMIAL  COEFFICIENTS  ARE  LISTED. 


E**** *************************** *************************************** 

DIMENSION  STR(/3)»  G( 7b) •  alPH(3)»  X(7b) 
oImENSIoN  eOlF(Iu) *  YI(75)»  SlfaHAT(A) 
uImlNSION  aNE* (75) »  YNek<7j) 
oImENSIuN  sT(1b) 

oA I A (ST ( I ) » 1=1 1 lo)  /.0b6> • oau*  O«u8o* .2bS» >255* .255» .5 tu* .52* .52* 

1  .7tl(i».7ua*.7G6i.8<*3».d43*.843»i.2bS*1.255rl.2b5  / 

REaD(5»**0£:o)  •« 
rNT=u 
iFLaG=(i 
i\FlaG=u 

U=A 

PRanT  2oOC 
t.0  TO  200 
aOu  KN l =1 

PRINT  2oOC 
oO  TO  3 JO 

ZOu  rEAd(  b#  17uO*eNu=4oO )  A1 »  A2»A3»  A4»  A5»  Ae>»  (ALPH(  I )  »  I=l»o) 
.aF(KNT.EQ.O)DTu=A1 
IF  ( Al.NE.  DTJ  )  GO  TO  500 
KN1=KNT*1 
300  TEMP=A1 

PRaNT  5b0o*  A4» Ab 

UH(=A2 

THEK=A3 

sTr(KNT)=S1 (J) 

G(KNT)=(ST(J)*aS)/A4 
j5u  COnT INUt 
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u=o+I 
i\£r=A6 
JTo=Al 
vjO  10  2u0 
hOu  1FLAG=1 
30c  uPIS=KNf 

UO  600  l=i*NPI:> 

X<i)  =  1.u-lOS(3Tk<1)) 
oCo  COnTlNUc. 

WRITE  (otlbOo)  ALPhtUHl *TlmP»ThCK 
wRllt  lofiyOU) 

1=1 

lTLMP=NPTb 

?0v  IF  11. oT. HEMP)  oC  TO  lOOo 

OOo  *RlTE  (o#£ioOO)  bIK(l)#STR(I+l) *STR(i+2)*X(l>*XlI+l)#X(I+2)*G(I>*G( 
ll+i)  #GU+2) 

1=1+3 
toG  TO  700 
10u0  COnl  INUl 

C******************************************************^*************** 
CAcL  OUTLR  (NPrS*X#G»Nc.*PTS*ANEW.YNbW) 

v,**.******************************************************************** 

Wl=0 

Uu3  NI=Nl+i 

C«***********»******«*****»*******«*%**«********************.*********** 
CAwL  CVrtto  (NltNE*PT8»XN£fc* YwEta»UEGR£>LOEF) 
l************ *************************************** ******************* 

aibHAT  tuD— 0« 

u  polynomial  calculation 
uO  1200  I=l»NE*Plb 
Y 1 1  i  )  =CuEF C 1 )  +COEF  (2 )  *XhEw  ( I ) 
lF(NI,Eu.<:)YIli)=YIll)+C0u:t3>*XN£*ll>**2 
iFiNl.L«.3)YIU)=Yi(l)+C0tF(H)*XN£W(I)**3 
iF(Nl.Eu.4)YI(i>=YKU+C0LF(5J*XNEW(l)**4 
5lGHATli'iI)=SlGri*T(Nl)  +  lYNLrfU>-Yl(l>)**2 
I2u0  CONTINUE. 

uF=uEWPIS«nI-1 
bloHAT  ( Ni  )  =SlGrlAT  t Ul )  /OF 
IF  (NI.toT.l)  Gu  TO  1300 
oO  TO  HOC 

IouO  iFlNl  .NE.  2)  oOTO  1350 

«RlTL ( N» 40u0 )  l HCK , UHT  t TEMP, uEKPTSr (COlF ( I ) r 1=1 r  3) 

.wFurE(W*40l0>  lXN£ta(I)»I=l*NEWPTS) 

WRITE (Nt 4010 J  (YN£*(I)»I=1»i,LI»PTS) 

1330  COuTINUt 

IF  (NI.lT.4)  Go  TO  U00 
I4o0  IF  UFLAG.lQ.D  GO  TO  1500 
GO  TO  100 
13o0  STyP 

IouO  FORMAT  t///2X»»F  =• »£lb,a»5x» *SIG» rll» »SQ  =» »E15.8»5X* »SIG« # lit *SQ 
1  =*»E15.8) 
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ji /uO  FOkHAT  1 3F 10  •  u  *  2F  lo  •  b*F  lu  •  u  *  tjX  *  3Ah  } 

ltiuO  FOurtAT  (//1X»o«h»  *  •  UROH  Hc.Ii>hT  0F»  •  »  *  *  TErtPEKATURE  ♦#F6»1 

i»'M  THiCKl.ESb*  »FA •  1  # 2n*  *  */ ) 

iyoO  FOnfaAT  11M  *8X»  'b^KEbS  LtWcLb' »21X*  ,l«o**COb(STKESS  t-EvEl_)  *  *25Xt3F)6 
l'b> 

<;OoO  fOmMAT  (1H  t 3(OF1G.4#5a)  ) 

^iu3  FOnriAT  ll«i  »3(f1u»4»25A)  ) 

2cuO  POtv*iAT  cin  »3(tH0.4*lsx)J 

c3j0  FO«MATUHl»lbX**OKI6INAL  um Ta‘»//) 

43uo  FOnrtAT ( JFb»l# jFl2*7) 

*♦013  POt\l*iAT(iuFo»4*/»lCFtj«4) 

H020  FUnHATlll) 
abuG  FOrtFiAT  t  dX»F  12«b*F  lb«5) 
ct\iu 
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SUokOUT INE  OUTlR  (NPTS,X*  Y  *NEwPTS*XNEW*  Y NE.M ) 

u* ********************** ♦******♦+*» ************************************ 

l 

l  OoTLlER  POINT  REJECTION  SubPROGRAm. 

c+ ********************* ********* **** ******** ******* ********* *********** 
l  ME  I  HOD  oASeD  ON  TnE  NAlR  CRITERION  (AN  EXTtNSlON  OF  TmE 
0  c-XThtME  STuu£nI12EU  uEvlATE  FROM  THE  ShMPLE  MEAN), 

C 

l  FOR  THE  SlT  OF  ObStRVAT IONS  hAVlNG  THE  LARGEST  VARIANCE 
l  tACH  ObsERvAT low  OF  THt  SET  .IS  TESTEO  INDIVIDUALLY  AS  A 
C  CANDIDATE  FOR  REJECTION  AS  AN  OUTLIER. 

C 

U ******************** ************************************************** 

ulMENSlON  X(l)»  Y  ( l) »  XNEh(l)*  YNEW(1>»  XA(25)r  XB(25)»  XC(25>»  YA 
1  ( 2‘j )  r  Yb(2b)  *  YC (2S) »  YMEAN(^5)  »  S2N(2s) •  REY(2S) 

TEblV  =  1./3** 

UO  100  1=1*25 
loo  is£Y(I)=0 

ImTc.MP=NPTS-2 

l=1 

1=1 

200  lOnTINUl 
J=l  +  1 
*=1+2 

IF  (I.GI.U1EMP)  faO  TO  oOO 
GoTO  25C 

IF  (X(J) .NL.XIR) )  GO  To  4uG 
IF  (X(I).NE.X(J))  GO  To  300 
£.bo  CuUTINUE 
XA (L)=X ( I ) 
xB(L)=XlJ) 
aC(L)=X(K) 

YA(L)=Y ( I ) 

YblL)=YU) 

YC IL) =Y (K ) 

*EY(L)=0 
L=L+1 
1=1+3 
GO  TO  20G 
30o  COnIINDl. 

1=1+1 
GO  TO  2u0 
*+0o  CONTlNUt 
1=1+2 
oO  TO  200 
aOu  COuTlNUe. 

L=l-1 

NP=L 

nPTS=(L*3)"1 

C  CALCULATE  means  G  VARIANCES 
oO  700  1=1»NF 


112 


AF  (YAU).hE.Yd(1>.0R.YA(1).kE.YC<1))  oO  TO  oOu 

YMuAN (IJ-Ya(X) 

b2N(i)=U.U 

rtRlTt  lo»£:oOu)  YMEaM  I )  t  S*:N(  1 ) 
bO  To  7u0 
uQo  LOnlINUe. 

YMtAN ( I ) - ( 1 A (1 ) +YO ( I ) + YC (1 > ) /3. 

b2iM(I)  =  ((Y«(I)“YMt.AN(I))**2+{YB(I)-YMEAN{I))**ii+{YCCl)-YMEAN(I))** 
1D/2. 

/Do  uOuTlNUe. 

C*  +  *  +  ***  +  ***************  *********,*  4-***  ***************  ******  ************** 

C  sOrtt  IN  OauEK  Of-  DECREASING  VARIANCE 

CAUL  SORT  <NP#S*;n*YMEAn»YA>Yi>»YC»XA(XB»XC) 

K.***** ****************************  ************************************* 

L  uAuCULA TE  SNU 

u=i 

C  SUi-i  VARIANCES 

bOO  RNf=U 

sUi»,V=0. 
oO  900  1=1 » NR 
If  (KEY(l).EQ.i)  00  TO  9uu 
if  (I.Lu.L)  60  10  90u 
sUhV=S0mV+S2N(1) 

KNi=KNT+1 
90u  OOnTINUl 

SNu=SGKI (suMVZ(RNT-I) ) 

C  TLaT 

TA=ABS(  (YAlL)-YrtE«N(L))/SiMU) 

TB=AbS< (YulL)-YtoEAN(L) )ZSNU) 

IC=AbS<  ( YC(U“YMEAN(L)  )ZSh(U) 

iF  (TA.cU. IB.OK.TA.Eui.TC.OR.lB.EQ.TC)  uO  To  1300 

TE=MAX(  TA*TR»  ID 

IF  (TE.lE. IES1 v)  oo  TO  13u0 

IF  (TE.uQ.IA)  bO  TO  lOUO 

If  (TE.c-G.TB)  bo  To  lluO 

IF  (TE.tQ. 1C)  bO  To  12uO 

bO  TO  1300 

xOoO  *RiTE  (o*2A00)  UrXA(L)  t  YACL)  rTErBNU 
.  aA  ID  =999* 
rEY(L)=1 
uO  10  1300 

iiuO  (/Rate  <o (Eaou)  l.xuil)  »yb<d  (Te^snu 
ABlL)=9y9. 
r\EY(L)=l 
bO  TO  1300 

lc:uO  *FuTE  (b»240o)  L»XC(L) »YC<D »TErsNU 
XCIU=999. 

REY (L)=l 

0  CHtCK  TO  StE  IF  ALL  VALUES 
13uO  IF  (L.tb.i/F)  GO  TO  140o 
l=u+1 
60  TO  800 
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FUI  IN  UE-w  AKrlAYb 
L=u 

i<EwHTS=0 

U=L+1  «  ,  ft 

it-  (XA(u)-y99.J  looO»l/Ou»louu 
FOiNT  ACCtpTtj 
NE»PTb=NE.*HTb+l 

aNc. «»  ( NE*PTb )  =XA  l  U I 

tNe.«lNe.AP<bi=YAU.) 

IF  lXBlt-)“999.)  Ibu0#ly00»  l'b«u 
uEaPTSSixEaPTS+l 
XN'cn  (Nt«PTb)—  Xo(Ul 
YNtw  (Nt-oP  fb)  =  Yt»U-) 

IF  (xClL)-999.)  20u0*2l0U»2uu0 

wE»PTb=n£*PTb+i 

Al')L.*i(Nt.nPTb)=XOUj 

YNi_a (NfcaPfb) =  YC (t) 

IF  lL.Oc.NF)  TO  2200 
u0  TO  lbOO 
Kt i uKN 

FO«i*AT  Uriu»12A»*bAMt  *ALUcd  **2Fb.O)  -»,F6.2.' 

FOkMT  I'GKEJtCl  POINT*  #  I^f  *  *  X  =*.Fd.Ar*»  Y  -*.Fb.2»'»T  - 
1 #  SuO  =  • » Fo • 2 ) 
cNu 
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bUuKOUT INt  CVKLb  (WirNPT# A» Y #OtOK£»COEF ) 

v^*********************^ 


********** 


v- 

u 

L 

C 


tSSitSSS  DiUKNl,^  the  EOUAUOU. 

OF  THE  bTATlSliOALuY  2eST  FlUlNb  POLYNOMIALS  OF  FIRbT » 

SEl0ND»  THIRL)  AND  FOURTH  ukuER. 


I.****************  *******  **** 

jIi*.eNSIuN  A(lo»ii>*  XU) 


*********** 


w=Nl+l 

M=l«+1 

anht=nht 

U(J  100  i=lrfl 
UO  luO  vl— 1  * M 
lOu  A(i»«J)=o. 

A(x»l)=wPT 
DO  400  n=i»npt 
jO  500  i=i*M 
uO  300  o=l»N 

!Po2=I+5-^ 

IF  (IPJE)  500»5GO*20u 
„Qo  A(i»U)=a(I.J)+a(KJ**IPJ2 
jOO  CONTINUE 

uO  ADO  1=2  »N 

40o  a(1»M)Sa(I»M)+a(K)**(I“1)*Y(kj 

uO  bOO  J=1»MPT 

L>0u  A(i»W)=AllfM)  +T  IU> 

UO  600  <J=1»M 
L)0  oOO  i=l*N 
oOO  A( x I » J) /hNPT 
iErcR=0 


/Ou 

oOo 


you 


lOuO 

IluO 

icuO 

i5G0 


I4u0 

1500 


M=N+1 

JO  130 J  l—i *  N 

IF  lA(l»l>>  8u0»7uu#a0u 

lEr\K=l 

uO  TO  1400 

TEi*iP=1  •  U/A 1 1 »  I ) 

IP1=1+1 

uO  900  j-IPI *  M 
A(1*J)=a<1»«J)*1eMP 
uO  1200  K=1*N 
IF  CI-K1  1000» i200» 1000 
DO  1100  J=iPl*iv' 

A  (A*  J)=a(K»J)-AIK*1)*AU»U> 

CONTINUE 

CONTINUE 

N=wl+1 

rt-N+1  , _ 

IF  T1ERH)  150u» 16uy» I5u0 
WRITE  to»3l0u> 
t.0  TO  2500 


115 


4t>uO 

1.700 


IcUlO 

19uO 


<1000 


cloO 


2200 

2.100 

c 

<1400 

2buO 

cbuO 

27u0 

<1600 

<1900 

30u0 

3100 


COuTINUl 
JO  1700  K=1»N 
COtF (K)=A(K»M) 

6UtoR2=0.G 

JO  1900  I=1»NPT 

YC=C0EF(1) 

JO  1600  K=2*N 

YC=YC+C0EF(K)*X(I>**(K-1) 

K=Y  1 1 J  — rc 

SUMK2=SUMRc+R*K 

SI6MA=S0RT(SUrth2/ANPT> 

SStRR=SUMK2 

SUMR2=Y<1> 

JO  2000  I=2»NPT 

6UmR2=SJMR2+Y(1) 

t>AKYl=SUMK2/NPT 

6UmK2=0*0 

DO  2100  I=1.NPT 

K-Y  (I)“6ARY1 

SUwR2=SUMR2+R*R 

SST0T=SUMR2 

SSKEG=SbT01 -SStKR 

JSR£0=SbRtb/Nl 

USt.RR=Sb£RK/ (NPT—  t  Nl+1 ) ) 

FRaT  0=USREb/DStKR 

JEbFT=Nl 

0EtoF8=NPT-(Nl+l) 

UEGR£=NPT-1 
c.Tb=SSERR/bSTOT 
iFlE1S.o£.1.0)tTS=1.0 
CORK=Si»RT  ( 1 .  O-ET  S ) 

write  (faooooi  >n 

wRiTE  (o»2b00) 

WRITE  (b»2700)  SSREG  . 0E6FT  . UbREG  »  FRAT 0 » CORR 
WRITE  Co. 2600 )  SS=RR  » DEGF 6 * USERR 
•RITE  (e.2900)  SST0T»0t6RE 
WRITE  (6.2400) 

MM=N1+1 

DO  2200  1=1.10 
K0im(I)=1-1 

WRITE  (o»2bOO)  (KOn(I) .COtF(l) .I=1»MM) 
RETURN 


FORMAT  (//»2X» 'CURVE  COEFFICIENTS') 

FORMAT  (//2X.2Hbl»ll.lH)«3X.tlS.7) 

F0RMAT(//»2X» 'SOURCE* »9X»  *S.S»  »9X» 'O.F'»9X-«  *M.S' i-9A» 'F  •  »12X.  »R»  > 
FORMAT  (//.EX.'DUE  TO’ »5(4X»£10.4) ) 

(//»2X. 'ABOUT' »3(4X»£i0.4) ) 

1//»2X.*T0TAE'»2(4X.E1Q.4>) 

(1H1.2X»»AN0VA  FOR  CURVE  OF  0RDtR'»I3) 

(//»4X. 'SINGULAR  MATRIX  * »/.4X» 'CURVE  FIT  IMPOSSIBLE*) 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
ENu 


j  J 


j 

i 
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SUbKOUT INE  SOKT (N» VAL»Xl#X2»X3»X4»X5»Xb»X7) 

(..Mm******************************************************************** 

0 

v,  SUbROUTlNE  FOR  SORTING  N  NUMbERS  IN  DESCENDING  ORDER 
L 

C*»* *********** ***************** **************************************** 

uIi'iENSION  vAL(i)»Xl(D  #X2(1)  »X3(1>  #X4<1) »Xb(l) *X6(1> »X7(1) 

W-IM-l 

UO  100  i-l»M 
Ul‘l 

UO  100  1I=L»N 

IFlVAL(l)  .GE.  tfAUil))  GoT'j  loO 
F  =  VAL(I) 

VAL(I)  =  VAL(Ii) 

VAu(II)  -  F 
F=Al(I) 

xiii)=xiai> 

Aiai)=F 

F=A2(I) 

x2a)=xaai) 

X2(II)=F 

F=X3(I> 

A3U)=X3(ii> 

X3(II)=F 

F=a4(I) 

A4(I)=X*+Ul) 

X4(1I)=F 

F=a6(I> 

X5lI)=Xb(ii) 

xsui>=f 

F=a6(I) 

X6(I)=Xq(Ia) 

Xb(lI)=F 

F=A7(I) 

X7lI)=X7(Il) 

A7UI)=F 
100  COwTINUE 
RETURN 
ENl> 


APPENDIX  E.  TOTAL  BOX  DATA  ANALYSIS  CODE 
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<.*«*♦***************#*************************♦*********************«** 

c 

c  TOTAL  BOX  UATA  ANALYSIS 

0 

c *********** *********** ******* ***************************************** 

c 

0 

C  MAIN  DRIVER  READS  THE  CUSHION  MATERIAL  DATA  RECORD  CONSISTING  OF 

c  temperature#  ukop  height#  thickness#  box  weights  #  g-value 

c  uOa  number  ANU  MATERIAL  TYPE  (IN  THIS  ORDER)  . 

L  INITIALIZES  ARRAYS  CALL  OUTER  (OUTLIER  SUBPROGRAM)#  AND  CALLS 

C  CVREG  (CURVILINEAR  REGRESSION  SUBPROGRAM). 

C 

c 

C  THt  F  -  STATISTIC  AND  THE  FIRST#  SECOND#  THIRD  AND 

C  FOURTH  oEGREE  POLYNOMIAL  COEFFICIENTS  ARE  LISTED. 

C 

C 

C***************************«*** ***********************************'*** 

c 

DIMENSION  WT(7s) #  G(7S) #  ALPH(3)»  X (75) #B0X(75) 

DIMENSION  COEF(IG)#  YI  (75) »  SI6HATI4) 

DIMENSION  XNEW(75) »  YNEW(7B) 

KNT=0 
IFuAG=0 
KFLAG=0 
GO  TO  200 
lOu  KNT=1 

uO  TO  300 

EOu  READ(  5#17uO»£nD=4uO)  Al »m2» A3# A4» A5» Afa# (ALPH(1)#I=1#3) 
IF(KNT.LQ.u)DTo=A2 
IF  (A2.NE*  DTJ  )  GO  TO  50u 
KN1=KNTU 
300  TEMP=A1 
Uhl=A2 
THCK=A3 
DTw=A2 
*T(KNT)  =  m4 
BOA(KNT)=Ab 
G(kNT)=A5 
uO  TO  200 
400  IFlAG=1 
SOU  wPTS=KNT 

DO  BOO  I=i»NPTS 
X(l)  =  wT(l) 
uOu  CONTINUE 

WRITE  (b»lt>OG)  ALPh»DHT»TeMP#ThCK 
WRITE  (O.I900) 
i  —  i 

itemp=npts 

70u  IF  (I.GT.ITEMP)  GO  TO  1000 

WRITE  (o#2100)  WTU)#  u(I)  »BOX(I) 
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aboKOUUNE  OUTLr  (NPTS*X*Y*NlWPTS*XnEW* YNE*) 

\.*4 ******************** *********************************** ************* 


V, 

L  OUTLIER  PoiNl  REJECTION  SUbPROGRAM 

v. 


<-  ME  I  HOD  uAStD  ON  THE  NAiR  CRITERION  (AN  EXTENSION  OF  THE 
L  I-XIREME  STuDENTIZE.0  DEVIATE  EROM  THE  SAMPLE  MEAN). 


w 

c* *********** **********+******** *************************************** 
DIMENSION  X(l)»  Y«l)*  aNEw Cl) *  YNEw(l) >  XA<25)»  XB<25)»  XC(25>»  YA 
ll2o)»  Yb(2a)»  YC (25) »  YMEANC25)*  S2N(25)»  KEY<25) 

TLbTV  =  1.734 
DO  100  1=1*25 
lOu  KEY ( I )=0 

NTthPSNPTS-2 
L=1 
1  =  1 

dOu  lOnTINUl 
o=i+l 
K=i+2 

IF  (l.GI.NTEMP)  GO  TO  bOO 
Eau  Continue 

xA (L) =X ( 1 ) 

AB(L)=X(J) 

aC(l)=xik) 

YA (L)=Y ( I ) 

YB(L)=Y(J) 

YC(L)=Y(K) 

KEr (L)=U 
L=L+1 
1=1+3 
GO  TO  20C 
aOo  CONTINUE 
1=1+1 
GO  TO  200 
40y  CONTINUE 
1=1+2 
GO  TO  200 
bOu  CONTINUE 
l=l-1 
NP=L 

NPTS=(L*3)-1 

C  CALCULATE  means  &  VARIANCES 
uO  700  1=1 *NP 

IF  (YA(1).NE.Yb(I>.0R.YA(I).NE.YC(I))  go  TO  60u 
YMlAN(I)=Ya(I) 

S2n ( I ) =0 . 0 

•RITE  (b*2300)  YMLAN(I) *S2NII> 

GO  TO  700 
oOO  CONTINUE 
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YMtAN(l)  =  (YA(I)+YB(I)+YCa))/3. 

S2M I  )  =  l ( Ym( I)-Ym£AN( I) ) **2+(YB(I)-YMEAN(I) )**2+(YC t I)-YMEAN(I) )** 

12)/2. 

700  CONTINUE 

C** ************************************ ******************************** 

c  sort  in  order  of  decreasing  variance 

CALL  SOKT  (NPf S2N*  YMEAN# YA*YB*YC*XA*XB*XC) 

c ********************************************************************** 

c  calculate  snu 

L=1 

C  SUM  VARIANCES 

BOO  KNT=0 

SUMV=0. 

DO  900  i=l#NP 
IF  (KEY(I).EQ.l)  Go  TO  900 
IF  U.EU.L)  60  TO  900 
SUMV=SUMV+S2Nll) 

KNT=KNT+1 
900  CONTINUE 

SNU=SQR T ( SUM V/ ( KNT-1 ) ) 

C  TEST 

T A=ABS ( ( Y A ( L ) - YMEAn ( L  > ) /SNU ) 

TB=A8S ( l YB (U -YMEAN  t  L) ) /SNU } 

TC= ABS ( ( YC (L ) -YMEAn ( L ) ) /SNU ) 

IF  (TA.tQ.TB.OR.TA.EQ.TC.OR.TB,EQ.TC>  GO  TO  1300 
TE=MAX(TA»TB*TC) 

IF  (TE.LE.TESTV)  Go  TO  1300 
IF  (TE.EQ.TA)  GO  TO  1000 
IF  (TE.EQ.TB)  oO  To  1100 
IF  (TE*LQ.TC)  GO  To  1200 
GO  TO  1300 

1000  wRiTE  (b*240Q)  L*Xa(L)*YA(L) *TE*SNU 
XAIL)=999. 

KEY(L>=1 
GO  TO  1300 

1100  NR1TE  (6*2400)  L.Xb(L) * YB(L) *TE*SNU 
XB(L>=999. 

KEY(L)=1 
GO  TO  1300 

1200  WRITE  (6*2400)  LtXC(L) * YC(L) #TE*SNU 
XC(L)=999* 

KEY(L)=1 

C  CHECK  TO  SEE  IF  ALL  VALUES 
1300  IF  (L.EQ.NP)  GO  TO  1400 
L=l+1 
GO  TO  BOO 

C  PUT  IN  NEW  ARRAYS 
1400  L=Q 

NEWPTS-0 
1500  L=L+1 

IF  CXA(L)-999.)  1600*1700*1600 
C  POINT  ACCEPTED 
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SUbROUTlNE  CVREG  (Nl*NPT*X»Y*DEGRE*COEF) 

C«**************************«*******«*************«******************** 

c 

C  CURVILINEAR  regression 

«.  curvilinear  regression  determines  the  equations 
C  OF  the  statistically  best  fitting  polynomials  of  FIRST* 

C  SECOND'  THIRD  AND  FOURTH  ORDER* 

C 

c«**»****************************************************************** 

DIMENSION  A<10*11>*  X(l>»  Y(l>*  COEF(l>*  KON(lO) 

N=N1+1 
M=N+1 
ANPT=NPT 
UO  100  I=1»N 
DO  100  J=1*M 
100  A( I*J)=0* 

A(l*l)=NPT 
DO  400  K=1*NPT 
UO  300  I=1*N 
DO  300  J=1*N 
1Pj2=1*J-2 

IF  (IPJ2)  300  000*200 
cOu  A(i*J)=A(I»J)+XtK)*»IPJ2 
300  CONTINUE 

DO  400  1=2  *N 

400  A(1*M)=A(I*M>+X(K>**<I-1>*Y(K> 

DO  500  J=1»NPT 
oOO  A(1*M)=A(1*M)+YIJ) 

DO  600  U=1*M 
DO  600  1=1 *N 
600  A(I*J)=A(I*J)/ANPT 
1EkR=0 
M=N+1 

DO  1300  I— 1»N 
IF  (A(IfD)  800*700*800 
700  IERR=1 

GO  TO  1400 

bOO  TEMP=1.0/A(I*I> 

IP1=I+1 

DO  900  J=IP1»M 
900  A(1*J)=A(I*J)*TEMP 
DO  1200  K=1*N 
.IF  (I— K)  1000*1200*1000 
1000  DO  1100  J=1P1*M 
1100  A(K*J)=A(K*J>-A(K»IJ*A(I*J) 

1200  CONTINUE 
1300  CONTINUE 
N=N1+1 
M=N+1 

1400  IF  (IERR)  1500*1600*1500 
1500  WRITE  (o*3100> 

60  TO  2300 
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1600  CONTINUE. 

UO  1700  K=1.N 
1700  COtF(K)=A(K*M) 

SUMR2-0 • 0 

UO  1900  I=1*NPT 

YC=C0EF(1) 

UO  1800  K=<:*N 

18uO  YC=YC+C0EF(K)*A(I>**(K-1) 

R=Y ( I)-YC 

1900  SUmR2=SUMR2+R*R 

SIGMA=SQRT  (SUMK2/ANPT) 

SSERR=SUMR2 

SUMR2=Y(1) 

00  2000  I=2*NPT 
2000  SUMR2=SUMR2*Y(i) 

8AKY1=SUMR2/NPT 
SUMR2=0.0 
DO  2100  I=1*NPT 
R=Y(I)-UARY1 
2100  SUmR2=SUMR2+R*K 
SST0T=SuMR2 
SSR£6=SST0T-SS£RR 
DSREG=SSREb/Nl 
OSERRrSSERR/ ( NPT- ( Nl+1 ) ) 

FRaT0=DSRE6/0SERR 

UEoFT=Nl 

0EGFB=NPT”(N1+1) 

0E6RE=NPT“1 
ETS=SSERR/SSTOT 
IF(ETS.6E.1.0)ETS=1.0 
CONR=SQRT ( 1 . O-ETS ) 

■RITE  (6*3000)  N1 
•RITE  (b* 2600) 

•RITE  (6.2700)  SSRE6.OE6FT.OSRE6.FRATO.CORK 
■RITE  (b.2b00)  SS&RR . DE6F6 . USERR 
■RITE  (b.2900)  SSTOT.OtGRE 
■RITE  (6.2400) 

MM=N1-H 

00  2200  1=1.10 
2200  K0N(I)=1-1 

■RITE  (O.2S00)  (KON( I ) .COLF ( I ) . 1=1* MM) 

2300  RETURN 
C 

2400  FORMAT  (//.2X. ’CURVE  COEFFICIENTS’ ) 

2500  FORMAT  (//2X»2H8(»Il.lH)»3X»El5.7) 

2600  F0RMAT(//.2X. ’SOURCE’ »9X» »S.S* »9X» *D.F» »9X» *M.S» »9X. »F » .  2X» »R* ) 

2700  FORMAT  (//»2X.’DUE  TO* »5(4X»E10.4) ) 

2800  FORMAT  (//»2X.’AB0UT’.3(4X.£10.4)) 

2900  FORMAT  (//»2X» ’TOTAL' »2(4X*E10.4) ) 

3000  FORMAT  (1H1»2X» * A^OVA  FOR  CURVE  OF  ORDER* »I3) 

3100  FORMAT  (//»4X» ’SINGULAR  MATRIX  • »/.4X. ’CURVE  FIT  IMPOSSIBLE*) 

ENu 


t  SUbKOUTINE  FOK  SORTING  N  NUMBERS  IN  DESCENDING  ORDER 
M=N-1 

DO  100  1=1»M 


L=I+1 

DO  100  II=U»N  • 

IF(VAHI)  .GE.  VAL(ID)  GOTO  100 

F  =  VALID 
VAL(I)  =  VALdl) 

VALCII)  =  F 

F=Xl(D 

XUDsXlUD 

XUID=F 

F=X2(I> 

X2U>=X2(I1> 

X2(ID=F 

F=X3(I> 

X3U)=X3(I1) 

X3tII)=F 

F=X4(I> 

X4lI>=X4UI) 

X4lII)=F 

F-X51D 

X5CI)=XS(ID 

XS(1I)=F 

F=Xb(I) 

X61I)=X6(IH 

X61II)=F 

F=X7(I) 

X7U)=X7(I1> 

X71II)=F 
100  CONTINUE 
RETURN 
ENU  .  ' 
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I********************************************************************** 

c 

u  INTERIOR  BOX  VALIDATION 

U 

U********************************************************************** 

c 

INTEGER  TYPEM 

COMMON  TYPEM<3) »NV .CONST.COcF (45) 

COMMON  TP*UH.TC.SS.GL»NVR.V(45) 

DIMENSION  X(500)*  Y(500)  .  YM<500)»  YPH50G)* 

1  YPU(500)*b<3)*YL(500)*YU(500) 

DATA  (TYPEMU) #1=1*3)  /  •MtNi'.'CEL  **«BOX  •/ 

DATA  (CuEF(I) .1=1*45)  /0. 0*0. 0.0. 0*7. 1865626*6*0. 0*-0.57911897* 

A  7*0. 0.-0.53894447.6*0. 0*30. U88255.  0. 0*6b. 557902*4*0 . 0  * 

b  4.1031526*5*0.0* 10. 219992*0. 0*-0. 22965288*2*0.0/ 

C 

KAhGE=1.0 

CONST  =  -14. 475703 

c 

C  PRINT  CUSHION  MATERIAL  TERMS  AND  COEFFICIENTS. 

PRINT  3999* CONST 

PRiNT  4000* (I*C0tF(I)*I=l*4S) 

C 

C  NUMBER  OF  COEFFICIENTS  FOR  THE  MODEL 
C 

NV  -45 
NA  =45 
C 

400  PRINT  1700 
C 

READ  (5rltoO0*ENl)=l500>  TC.UH.TP.NPTS* lb(K) *K=1 *3) 

NPTA=NPTS+1 

REaD(5*300u)  (XILL) *LL=2.NPTA) 

READ (5* 3000)  (YILL),LL=2*NPTA) 

500  PRINT  2500*  TYPEM(l) »TYPEM(2) »TYP£M(3) *  DH*TC*TP 
X(1)=NPTS 
Y(1)=NPTS 

CAUL  C0NFIU(X*V*NPTS»Yu.YU»6) 

UELTA=0.05 
X(2)=0.088 
I  =c 

oOO  CONTINUE 

IFlX(I)  .GT.  2.5)  GOTO  700 
C 

C  SECOND  ORDtR  POLYNOMIAL 
C 

YU)=B(l)+b(2)*U.Q  -  COStXlI))) 

1=1+1 

X(l)  =X(I-1)+UELTA 
GOTO  500 
700  CONTINUE 
X(1)=I— 2 
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m>=w 

0O0  COnTINUl 

R=i-1 

uO  900  l=2>K 

Y(L)=Y(L)+b<3)*<l»0  -  C0S(X{L)))**2 
c 

L  oEnERAL  CUSHION  MATERIAL  MOuEL 

c 

SS  =  X(L) 

i,  **************  ********  ***********  ********  ******************  *********** 

CAUL  MODEL 

^,************t ********************************************************* 

YKCL)  =uL 
900  COnUNUL 
C 

C  CAlCULAIE  prediction  limits 

L 

NPTA=X(1) 

£**************»***************************************************»*** 
CALL  PREDIC ( X  #  Y  *  NPTA  t Yl  »  YU  *  b  » V PL  > YPU ) 

i,********************************************************************** 

C  FInO  MINIMUM  iucc  g-level 

c 

YMiN=1000»0 
DO  1200  1=2 »K 

IF  (YMlN.LL.YCi))  GO  TO  1200 
YMiNsYd) 

M=1 

1200  COnTINUl 

c 

C  DETERMINE  VAHu  model  range  from  bounded  IUCC 
C  and  PREDICTION  limits. 

c 

*MIN=X(H) 

xl=amin-range 
xu=xmin+range 
wRiTE  (of 2500) 

WRITE  (of 2b00) 

UO  140C  1=2  *K 

N=5H 

i8=3H 

IF(X(I) .GE.XL.AND*X(I).LE.XU)IB=3H  ** 

'  IF l YPL( I) .faT.YM(I) .OR. YPU( I) .LT.YM(I) )N=3H  * 

IF  (YPL(I).LE.O.O)  GO  TO  1300 

WRITE  (0*2 700)  X(I)»Y(I) *YPL(I> *I0»YM(I) *N»YPU(I) 
bO  TO  1400 

C  NEGATIVu  G- VALUES  SET  TO  -  . 

13u0  WRiTE  (b*2o0U)  X(l)cY(I)*iU»YM(I) »N» YPU(I) 

I4u0  COnTINUl 
l 

C  PRxNT  5000»B(l)»B(2)»e;5) 


i 

| 

I 
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3iluO  POnMATtaA*  'BID  =**F12.7*3X*'B(2>=  •  *Fl2.7*3X»  »B<3'=  »»F12.7) 
c  NEaT  CASE 
00)0  4uU 

«, 

L  ENu  OF  oOa 
AboO  wRiTE  (6*2900) 

CAuL  EAaT 
C 
L 

louO  FORMAT (3Fb.2*  a3*3F12.7) 

17o0  FORMAT (1H1) 

tluO  FORMAT  l/7/*22X*»0«»Flb.8> 

22uO  FORMAT  (18X*Ib*FAb.8*5X*lbA4) 

23jO  FORMAT  (1H1*13X*3A4*4X*F4.1* •  IN.  D.ri.  • *F7.1* •  IN.  ThICK* *F8.1» • 
ITEmPERA  TURL  * ) 

2boO  FORMAT  (///16X* ’STATIC  SThtSb’ * 17X* ‘ACCELERATION  (G)‘) 

2t>00  FORMAT  (21X* ‘Pbl • * 13X* • IDCC»  *4X*  * LOwER-P • *  9X  * • MODEL • *  bX  *  • UPPER-P • ) 
2700  F0rMATCI8X*F7.4*9A*F7.2*4X*P7.2*4X* A3*F7.2*A3*3X»F7.2) 

2dOu  FORMAT (I8X*F7.4*9X*F7 *2 *7X»  *—  ** *5X* A3»F7.2* A3*3X*F7.c> 

2900  FORMAT  (1H1*‘  twD  OF  JOb») 

30o0  FORMAT tlOFb.4) 

3999  FORMAT ( AH1 *  5X  » • CONST • #  F 15 . 7 ) 

4000  format t 45 tax* ij»fis»8*/)) 

ENu 
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bUoROUT iNE  CONF ID  l  X  »  Y  »  NPTS  * 1 L  *  YU * B) 

c*»»******************************************************************* 

C 

L  COMPUTE  THE  PREDICT ION  LIMITS  OF  THE  CURVES 

C 

C*******************«************************************************** 

COMMON  TYPtM 1 3) *  NV  *  CONST  *  COEF (45) 

COMMON  TP*DH*TC*SS*6L*NVR* V(45) 

DIMENSION  X(l)*  tU>»  YPUll).  YPL(l)  »  till  >  *YL(l)  *YU(1) 

DIMENSION  aARIsOO'3) *  YAR(5oO>*  C(3)*  A(3*3>*  XIN<3*3)*  E(3) 

TAh=  1.734 
OX=0.0 
TS=0.0 
NPlS=X(l) 

DO  300  1=1*3 
Cll)=0.0 

tuco.u 
UO  100  t-=l*3 
XIn(I»L)=0.0 
A(l*L)=0.u 
100  COnTINUL 

jo  200  U=1*NPT5 
XAhU*I)=0.0 
YAKlU)=0.0 
cOO  CONTINUE 
30u  CONTINUE. 

F=u«0 

S=0«0 

SS«=0.0 

UO  400  1=1* NPTb 
0=1+1 

XAK(I*1)=1.0 
XAKlI*2)  =  XU) 
aAKII.3)=XARII*2)**2 
YAr(I)=YU) 

400  CONTINUE. 

DO  600  1=1*3 

UO  500  0=1*NPTS 

CU)=Cll)+XARU*I)*YARlJ) 

500  CONTINUE 

BX=  aX+LlI)*BU) 
bOO  CONTINUE 

UO  700  U=1*NPT5 
YS=YS+YARU)«+2 
700  CONTINUE 

SS«=(YS-BX)/(Xll)-3.0) 

S=5QRT(SSGi) 

DO  800  U=1»NPTS 

All»l)=A(l»l)+XARU*l)*XAKU*l) 

All*2)=A(l*2)+XAR( J*l)*xARlJf 2) 

A(1*3)=AU*3)+XARU*1)*XAkU*3) 

AU*l)=A(2*l)+XAR(0*2)*XARU*l) 


A(2»2)-A  (2#2>  +AAR  l«j*2)*XAK(  j»2) 

A(^»3)=ft(2f3)+/.AKlu»2)*XAKlOf  3) 

AU»l)=M(3r l)+XAR(J,3)*XAK(o*l) 
a  ( it  2  J  =A  ( 3*2-'  ♦AARlU*  3)  *XAK  (0*2) 

A(0f3)=A(3».5)+AAR(J*3)*XAK{j#3) 

COuTIWUt 

u=Mil»l)*{Ml2»^)*A(3»3)-A(3»2)*A(2»3) )+A(l»2)*(Al3»l)*A(2»3)-A(2*l 
1)*a(3»3)  ‘♦A'lO>*lA<2»i)*A(3*2>-A(3»l.>*A(2»2>> 

Xli't(l#i)  =  (A.‘2*2)*A(3»3>-Al3»2)*A(2»3)  )/D 
a1h(1»2)  =  (A(3»2)*A(1,3)-AU,2)*A(3*3)  )/D 
AlNll»3)=(A(l»2)*A(2,3)-A(2f2)*A(l»3) )/U 
XlU(2'l)=(Ai3'i)*A(2'3)-At2fl>*A(3’3))/D 
Xlw(2#2)  =  <A(l»i.)*A(3»3)-A{3»i>*All»3>)/D 
xIn(2*3’~*a<2»1) *A(l#3)-All»l>*A(2*3) )/D 
xIN(3*1)-Ia(2*H*A(3»2)-A(3»1)*A(2»2)  )/D 
XIwC3»2)=(a(3»1)*A{1»2)— Atl»i>*A(3»2) )/D 
Xlu(3»3)=(A(l*l>*A(2*2)-Al2*i)*A(l»2) )/U  . 

RE  I URN 


EN TRY  PhEUIC (  X  r  Y  » t»PTS » YL *  YU , b » YPL r  YPU > 


uO  900  i=l*NPTb 
XAK(I»1)=1. 

J=i+1 

XAK( 1»2)  =  1.0-COS(X(J)) 
XAt\(I»3)=XAR(I»2):!'*2 
90o  CONTINUE. 

YPu(l)=NPTS 
YPu(l)=NPTb 
uO  1400  J=l»NPTb 
oO  1100  I=l»3 
UO  1000  K=l*3 

tU)=£U>+XAR(J»K>*XIN{K»I>' 
lOuO  CONTINUE 
1100  CONTlNUt 

UO  1200  1=1 t 3 
F=F+E(I)*XAR(J#I> 

1200  CONTINUE. 

F2=1.0+F 

F2=SQRTIF2) 

F=bURT(F) 

N=«J+1 

YU(N)=Y(N)+TAH*S*F 

YLlN)=Y(N)-TAH*S*F 

YPo(N)=Y(N)+TAH*S«-F2 

YPL(N)=Y(N)-TAh*S*F2 

F=0.0 

UO  1300  1=1 r 3 
£<i)=0.0 
1300  CONTINUE. 

14u0  CONTINUE 
RETURN 
ENU 
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SboROUTlNt  MOOtt- 

C *************** *********************** 

lNrtGEK  TYrEM  t**********************************^,^^ 

COMMON  1 YPtM (0)  *NV  r CONST » COcF (45) 

COciMON  1 P  #  OH » 1 C  #  SS  t  GL  *  N VR  t  0  ( t*5 ) 


c 


ML  =  l.O-CoS(SS) 

A La  -  At  *  AL 
SRbH  =  bORT (  OH  ) 
fR  =  <Tf'  +  46U.O)/iOG. 
(Hi  =  Tk  *  TR 
IKo  =  Tn  *  TR2 
Tk4  =  TR3  *  TR 
rCOH  =  1C  **  (-0.5) 
TCTH  =  1C  **  (-1.5) 


VIU1J 

V  ( 02) 
V(uO> 

V  ( 04) 

V  ( 05) 
1/(06) 

V  ( 07) 
v(0  0) 


TR 

TR 

TR 

TR 

TR 

TR 

TR 

TR 


*  TCOH 

*  TCOH 

*  TCqH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCOH 

*  TCOH 


SruH 

srdh 

sroh 

5K0H 


V(o‘4»  =  TR 

*  TCuH 

*  SRbH 

V(10)  “  TRa 

♦  TCOH 

V(ll)  --  TRa 

*  TCOH 

V<12)  =  TRa 

*  TCOH 

v (13)  =  TRa 

*  TCTH 

*  SRDH 

V(i4)  -  TRa 

♦  TCTH 

*  SRDH 

V(lb)  =  TR2 

*  TCTH 

*  SRbH 

Y(16)  =  TRa 

*  TCOH 

*  oRbH 

V((17)  =  TRa 

*  TCoH 

*  SROH 

V(ia)  =  TRa 

♦  TCOH 

*  SRbH 

V(19)  =  TRO 

*  TCOH 

V(aO)  =  TRO 

*  TCOH 

V(al)  =  TRO 

*  TCOH 

V(aa)  =  TRo 

♦  TCTH 

*  sroh 

V(aO)  =  TRO 

*  TCTH 

*  srdh 

V(a4)  =  TRO 

*  TCTH 

*  srdh 

V(a5)  =  TRO 

*  TCQH 

*  srdh 

V(t6)  -  TRO 

*  TCQH 

*  srdh 

V(a7)  =  TRO 

*  TCOH 

*  srdh 

V(a6)  =  TR 

*  TCTH 

V(29)  =  TR 

*  TCTH 

y'(00)  =  TR 

*  TCTH 

V(01)  =  TR2 

*  TCTH 

V (02)  =  TRa 

*  TCfH 

V(00)  =  TRa 

*  TCTH 

V (04)  -  TRo 

*  TCTH 

*  AL 

*  AL* 

*  AL 

*  AL2 

*  AL 

*  ALa 

*  AL 

*  ALa 

*  AL 

*  AL2 

*  AL 

*  ALa 


*  AL 

*  ALa 

*  AL 

*  ALa 

*  AL 

*  ALa 

*  AL 

■  ALa 

*  AL 

*  ALa 


*  AL 

*  AL«d 


V<3b>  -  TK3  *  TC-Th 

V  <  Jb)  =  TK.S  *  TCTH 

V(i7)  =  TR 
W(i8)  =  TH  *  Au 
V(39)  s  TK  *  AL2 
VUO)  =  TK*: 

V<*+1)  =  TKi.  *  AL 

VU2)  =  TR2  *  AL2 

V(43>  =  TK3 

V  (44)  =  TKb  *  At. 

Vl45)  =  TRd  *  AL2 

L 

t-  COMPUTt  DYNAMIC  CUSHIONING  FUNCTION 

c 

GL  =  CONST 
UO  100  U=1»NV 

lOu  oL  =GL  +COtF  ( J)  *  V  (U) 

RE (URN 
£NU 
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L  **********  **********4.*  ************************************************* 

V* 

l  TuIaL  box  VAlIoATION 

l 

C***********»*****************************************4‘***************4‘* 

COMMON  A <2. 3)  .COSH  (2)  .RON (2)  .Nl.NPT.  ILYA < 72) 
COMMONSt>K£o.Dbe.KR»T»ULlbl)  *UU< 51)  .ALU19)  .BARX 
COmMON  Y(lt>0)  .X(2'15u) » IEKR.  TC.DH.TP 
Nl  =  1 
T  =  1.704 
Rlm0(5»99)  ILYA 

10  RtAU(5>looO»£NO  =  9y9)  Dti»TC»TP»NPT 
uO  Oil  I  =  i*NPT 
KEmU(5»170u)  XU.l).Y(i) 
oO  tONflNUt 

*RlTE(b»89)  ILYA 
*RiT  t.  ( b»  1800 )  Ori.TC.TP 
uO  40  J  =  1  *  NP I 
40  *Ri  TE(6»79)  X(1»J>.Y<J) 

lalL  LRlG 

IF ( IERR  .Eo.  1)  GOTO  lo 
WrtIT£(o.o9)  Ilia 
mRiTE (6. 16uO )  Ofl.TC.TP 
Call  LIMIT 
uOTO  10 

yyy  wRi1E(b»19oO) 

call  EXIT 

c 

79  F0rMAT(1Ua.F10.4»H0.4) 

09  FORMAT  (1H1. 12X.  72Ai) 

99  FORMAT (72A1) 

iouO  FORMAT  (3Ft>.  2 » 10) 

1700  FOrMAT (32X»F8«4»5a»F8»4) 

1800  F0rMAT(13X.  <+X.F4.i.»  IN.  D.H.  *»F7.1.*  IN.  THICK‘» 

A  F8-1*'  TtMPERATURE*.//) 

1900  FORMATdHl.lOX. 'ENU  OF  JOB') 

ClMO 


c***+ 

c 

c 

c 

<-***♦. 


bUoHOUTlNL  LHtu 

********4.*,*****^ 


un-ftAh  KE&RESSioN 


*******  ******* 


<-0i.ii>i0NSbKEb*UbtRR*T,uL'5l)  ,,,.  l<ii  )  ai  77^.-,, 

2“;  ra50)**(2,15u).iEK^T cf^:?pai9,'OARX 

M=N+1 

aUPT=NPT 
00  IOC  i=l#N 
oO  100  J=1»M 
'  «(a*J)=u, 

A(x#1)=npT 
OO  400  A=1»NPT 
OO  400  1=1»N 
oO  3u0  U=1»N 
IPo2=I+u—2 

IF  (IP J2)  o0u»J0o»ci0u 
Ali»J)=A(I»0)+A(l»K)**ipJi> 

COnTlNUt 

UA«A  =  A  ( 1  *  M )  /Ai'JpT 
OO  400  i=<i#N 

?nlZ%l  ~  Aa',n,+x(A.K)**a-i)*y(K) 
oO  bOO  J=1»NPT  ; 

OO  600  o=l »M 
OO  oOO  1=1 » N 
A ( x  »  J } =/* ( I r  J ) / ANp  T 

IEaR=0 

M=W+1 

00  1300  1=1  *14 

IF  (A(ItU)  8oo*7uu*80o 
IEaK=1 
00  TO  1400 
rENtP=l,u/A(l,j; 

IPl=l+l 

UO  900  u=IP1,m 
A ( i > J) =a ( I r J) * JEmP 
00  1200  K=1»N 
IF  (I-K)  l000»12u0,lo00 

00  1100  J=lPl#rt 

fn^Tk,=:A(K'vJ,"AlK’1,*A<I''I> 

COnTINUc. 

CONTINUt. 

N=n1+1 

M=n+1 

IF  (IEKR)  1500*lb00»150fl 
*RlTE  (o»3l00) 

O0  TO  2300 
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XbuO  COwl INUt 

UO  1700  K=l*N 
i7v»0  COc.F(K)=A(K*M) 

SUhR2=0.0 

uO  1900  I=1*NPT 

YC=COEFtl) 

1)0  1B00  K=2*N 

itiUO  YC  =  YC+COc.F(K)*X(l#I)**{K-l) 

R=Y<1)-TC 

i900  SUrtR2=SUMK2+R*R 

SIdrtA=SttRT(S'JMK2/ANPT) 

SStRK=SUMK2 
SUmR2=Y(1) 
uO  2000  I=2*NP  f 
2O0O  sUhR2=S0MR2+Y U ) 
oAKYi=SUMKc/NPr 
5Uwr2=0.0 
UO  210U  I=1*NPI 
K=r (I)-bARYl 
2I0O  aUrtK2=SUMRt+R*K 
SSTOT=SuMKc 
sShEG=SsT01-Sse.RR 
uSr£G=SsREo/N1 
uStRR=SbERR/(NPT-<Nl+l) ) 

FRATO=DSRto/DStRR 

UEoFT=Ni 

uEoF6=NPT-lNl+l) 

0EGRE=nPT-1 
tTs=SSERR/sSTOi 
IF(ETS.GE.1.0)cTS=1.0 
COhR=SGKT  ( 1  •  0“C.TS ) 

•RITE  <b*3G00)  W1 
•  Rift  (b*2bOu) 

•RITE  (o»270u)  SSREG»QtGFT *OSREG*FRATO»CORR 
•RITE  (o*2b0u)  SSERR*Ot.GFd*uSERR 
•RITE  (b#2900)  SSToT*OtGR£ 

•RITE  (b*2400> 

MM-N1+1 
uO  2200  1=1*2 
2cOJ  K0nII)=1-1 

•RITE  (o*250u)  TK0NU)»C0EF{I)*I=1*MM) 

2300  .RETURN 
C 

2400  FORMAT  T//.2X* 'CURVE  COEFFICIENTS*) 

2500  FORMAT  (//2X*2Hd(*U*lH)*3X*£l5.7) 

2oU0  FORMAT (//*2X*  *SUURCE*  *9X*  *S. S’ *9X*  *D.F*  r9X*  »M.S*  *9X*  *F • , 12X* *R* ) 
2700  FORMAT  T//*2X**UUE  TO* * 5<4X.tlQ.4> ) 

2000  FORMAT  (//*2X» • ABOUT* *3t4X*El0.4) ) 

2900  FORMAT  (//*2X**TOTAL**2(4X*E10.4)) 

oOOO  FORMAT  UH1»2X»'AN0VA  FOR  CURVE  OF  ORDER’ *13) 

3100  FORMAT  ( // *  4X  *  *  SINGULAR  MATRIX  » */*4X* ’CURVE  FIT  IMPOSSIBLE*)  . 
ENu 
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bUaROUT INc.  Llriil 

c** ******************** 


****************  ******** 


CUl4FIuE.NCt  LlMl  1 S 


L*************************^^^+^^^^^  t 

COrihOI-j  M(2,3).i;0EF(2),K0N(2),Nl,MPT.ILrAC7ii)*********************' 

LOrtKONSbRtb , ObtKR *  T » UL { 51 >  » uL  « 51 >  » «L  < 1 1 9) ,  bARX 
LOiv,hON  f 1150)  »Al2»lbo)»I£KK,7C>DH#TP 
ulmtNSIuN  LM<5)  *  YR(5i) 

uArA^(CM(l).I=i,5J  /-6,  ,J4b74o7»  0«0*o.  07231763  »**0.  Qu9l  702196*0.0  / 

bSc.KR=SuRl  (DStKtO 
5UiiK=SSK£G/  (CocF  (2)  **2) 
l/ElTA  =  0.5 

X(i*l)=55.u 
AO  1=1+1 

X ( 1 » 1 ) =X ( 1 * I-i j +uElTA 
1F(a(1*I)  »l£ •  62.0)  oOTO  10 
NX  =  l-l 
M2=SURT(DH) 

1 H= ( TP+460 .  0 )  /l  Oil  •  u 
TM2=  TH+TH 
THo=TH2*  Th 
uo  20  1=  A, NX 

t ( 1 ) =COtF ( i ) +COtF < 2 ) *X ( l » i ) 

TER  =  H2+X ( 1  *  1 J 

XYM<icM7j:{^:^(2,*TH"TE^^^TH^T^C«<->-TH3.TCK 

;o  continue 

L»0  3u  i  =  a»  NX 

SfcS'rt^UwilwRl07  WT* 1 1 1 X 1 1 "  >  -«NX»  /SU»K ) . 

ULU)=YU)+SH£rti 
>0  ULlI>=YU)-SHEKI 

#RlTE(6*5Q) 

UO  40  1  =  1*NX 

0  C0NTlNUt51>  XU,1,fY(l>'Dt-(^^M(I,,UL(I) 

RETURN 

0  AF0KMA5^PP;S^i*  *  16X*4HI£BL*  16X*  5HLQWtR* 15X*5hM0DEL* 15x* 

1  .  FOkMAT(5UuX*F10.4)/) 

ENu 
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C****************V*********** ***************** **************** 

C 

C  total  BOX  VS.  iWTtKlOR  BOX 

G 

C ****************************** ******************************** 

INTEGER  TYPEM 

COMMON  1 YPLM(i) # NV # CONST *COEF( 45) 

COMMON  TP#UH#TC#SS#GL»NVR»V(45) 

J I MENS I oN  2(6) #CM(5) »Xll0u> t Y C10O) * YM(100) #YC(100) 
l)A1A  (TYPtM(I) »I=1#3>  /  »MINI*»»CEL  #.*BOX  •/ 

O/' ]  A  (COEF(I)  » 1=1  #45)  /0.0#G.0»0.0»7. 1865626 * 8*0 •  0 » -0 •  57911897 » 

A  7*0«U#— 0. 838^4447 »6*U«0» 30 .038255#  0.0. 6b.5579u2*4*0.0# 

b  4. 1031526# 3*0. 0# -10. 219992 »u.0#-0.22965288»2*0.0/ 

OATa  <CM(I)»1=X»5)  /-6. 3487467 » 0.0» 0.07231763 »“0. 0091702196# 0.0  / 
UA.-A  (Z(I)»I=1»6)  /0 • 088# 0.<:55r 0*520 #0.708# 0.843# 1.25b  / 

COnST  =  -14.^757oC 
C 

C  PkINT  CUSHION  material  terms  and  coefficients. 

PRINT  3999#C0N3T 

PRINT  4000#  (I.COEFd)  »i=1.43) 

L 

C 

C  NUMBER  OF  COEFFICIENTS  FOR  THE  MODEL 
C 

NV  =45 
C 

40o  PRINT  1700 
C 

DELTA  =0*3 
A(l)  =  35.0 
i=.l 

oOo  CONTINUE 

IFU(I)  .fit.  62.0)  GOTO  700 
L 

1=1+1 

XU)  =X(I-1)+ULLTA 
oOIO  600 
700  CONTINUE 
K=l-1 

800  CONTINUE 

HEAD  (S»1o00#LNwI=1500)  TC»Dh»TP 
H2=StfRTlOH) 

TH=(TP+460.0)/10l>.0 
TH«;=  TH*TH 
TH3=TH2*TH 
uO  900  L=I#K 
TEk  =  H2*XIL) 

C  tXTERIOK  BOX  MuDEL 

YM(L)  =  CM  ( 1 )  +4.4  ( 2 )  *TH*  TER+CM  ( 3  )  *TH2*TlR+CM  ( 4 )  * 

L 

L  GENERAL  CUSHION  MATERIAL  MOotL 


149 


c 

SS  =  X(u>/i72.S 

C*+*^******************* *********** **************************** 

CALL  MOuEL 

C**************************************************»*********** 

Y(L)  =oL 
90o  CONTINUE. 

c 

wRlTL<to#2J00)  l  fYREM(I)  #1=1.3)  #DH»TC#TP 
•RiTE(6»2b00> 

UO  1040  1=1 #K 
C  CUSHION 

YCll>  =  YM( I)  -  Y(l) 

«KxTE(6#51)  XU)»Ym(1)»YU)*YCC1) 

1040  CONTINUE 
OOTO  800 

ISoO  *Ri.TE(b»  1700) 

CAl-L  EXIT 

2300  FORMAT  11H1#13a#3A4#4X#F4.U*  IN.  U.H.  *#F7.1#*  IN.  THICK* #F8.1» * 
lTErtPERATURt**//) 

itjbO  FORMAT (15X#*Uua  XC.1GHT  *»2A»*  TOTAL  bOX* »5a» *  INTERIOR  bOX*# 

1  5X» ‘CUSHION*#/) 

si  FORMATt  ,.0X»Fi3.4»3x»Fl2.4#sx#Fl2.4»5X#F9.4#/> 
lbOO  FOrMATUFS.2) 

1700  FORMAT 11H1) 

'999  FOkMAT ( 1HI # 5X # • CONST «  » F 15. 7 ) 

4000  FOkMAT l 45 l LX# I3*F15»8*/)) 
tNu 


r.  c.  <■■ 


U 


SUbKOUTINt  rtOOtu 
c*********************** 

**************************** 


************** 

********************************** 


u***************** 

^SsiSESi^?" 


COMMON 


v. 


o 


!\L  =  l.u-CoSISS) 

*Lc  =  AL  *  AL 

SRuH  =  bQKTl  0«  >ft  . 

TR  -(TP  ♦  46u»U)/lOO« 

TR2  =  TK  *  TR 
TRo  =  TR  *  7K2 
TRA  =  TK3  *  TR 
TCOH  s  TC  ** 

....  .  —  rf'  AA  (  ^  1  •  t>  ) 


VlOl)  =  TK 
V102)  =  TR 
Vlu3)  =  TK 
vlv>4)  =  TR 
VIG6)  =  TR 
wlub)  =  TR 
V 107)  -  TR 
Vlub)  =  TR 
Vlu9)  =  TR 
VUO)  -  TR2 

Vlll)  =  TKii 
VI 12)  =  TR2 
V  U3)  =  TR2 
vim)  =  TR2 
VUb)  =  TRa 
Vllfa)  =  TRa 
V ll7)  =  TR2 
VUB)  =  TR* 
VU9)  =  TK3 
V 120)  =  TRJ 
Vl2l)  -  TRo 
Vta2)  =  TR3 
Vla3)  =  TR3 
vim)  -  TRo 
Vla5)  =  TR3 
Vlab>  -  TR3 
V 127)  =  TR3 
Vta6)  =  TR 
Vla9)  -  TR 
VUO)  =  TR 
V 131)  =  TR2 
VI 32)  =  TR2 


*  7C0H 

*  TCOH 

*  TCOH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCOH 

*  TCOH 

*  TCOH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCTH 

*  TCTH 


*  SRLiH 

*  SRDH 

*  SrlUH 

*  r»KL>H 

*  SKOH 

*  S>KDH 


*  SkQH 

*  bKuH 

*  SKQH 

*  arlUH 

*  SROH 

*  SROH 


*  SROH 

*  SROH 

*  SROH 

*  SROH 

*  SROH 

*  SROH 


* 

AL 

* 

AL2 

* 

AL 

* 

AL* 

* 

AL 

♦ 

AL2 

* 

AL 

* 

ALa 

* 

AL 

* 

AL2 

*  AL 

*  AL2 

*  AL 

*  AL2 

*  AL 

*  AL2 

*  AL 

*  AL2 

*  AL 

*  AL2 

*  AL 


r  o  r> 
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vOi)  =  TR*  *  TCTH  *  AL* 

V(^4)  =  T Ri  *  TCTH 

V(jb)  =  T*0  *  TCTH  *  AL 

v(J6>  =  TKb  *  TCTH  *  AL* 

V(o7)  =  Tk 

V(btt)  =  TR  *  AL 
=  TR  *  Al2 
W(40)  -  Tk* 

V(4l)  =  TR*  *  AL 

V(42)  =  Tk*  *  ml2 

V(4i)  =  TtO 

*(44)  =  HO  *  ML 

VU5>  =  TfO  *  AL2 

CoMHUTL  UtNAMlC  cushioning  function 

GL  =  CONST 
uO  ICO  J— 1»NV 
100  GL  =GL  +C0EF{J)  *  V<J> 

KE I  URN 
fc.NL 


